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16.  Abstract 


Diffuse  multipath  is  a random  phenomenon  that  arises  out  of  the  diffraction  of  waves  from  rough 
surfaces.  This  form  of  multipath  can,  in  turn,  affect  the  angle  estimation  performance  of  monopulse 
antenna  systems.  In  order  to  evaluate  its  effect,  it  is  necessary  to  obtain  the  channel  spread  function 
that  results.  This  function  is  basically  the  wave  number  spectrum  of  the  resulting  random  field 
generated  by  the  scattering.  To  do  this,  use  is  made  of  the  Kodis-Barrick  scattering  cross  section 
and  Wagner’s  shadowing  function.  The  evaluation  is  performed  for  the  specific  geometry  of  a 
spherical  earth.  Evaluation  of  the  system  performance  degradation  is  also  made  for  various 
system  configurations. 
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Section  1 


SUMMARY  AND  RESULTS 

When  a signal  propagates  from  a source  to  a receiver  over  the  earth’s 
surface,  the  field  incident  on  the  receiving  aperture  is  composed  of  both  the 
free  space  component  and  those  components  arising  from  reflections  and 
diffractions  from  the  interface  between  free  space  and  the  earth.  The  reflected 
and  diffracted  component  is  commonly  broken  into  two  separate  parts:  the 
specular  multipath  return,  and  the  diffuse  multipath  return.  The  specular 
component  is  usually  well  defined  in  terms  of  its  amplitude,  phase,  and  incident 
direction.  It  is  a plane  wave  and,  as  such,  is  uniquely  specified  by  its  direction. 
The  diffuse  component,  however,  arises  out  of  the  random  nature  of  the 
scattering  surface  and  as  such  is  nondeterministic.  It  does  not  come  from  a 
given  direction  but  from  a continuum.  It  is  this  form  of  multipath  that  is 
studied  in  this  report. 

1.  1 Problem  Definition 

The  earth’s  surface  has  a height  distribution  that  is  random  in  nature 
so  that  when  one  tries  to  solve  the  propagation  of  electromagnetic  waves  over 
such  surfaces,  the  solution  must  take  this  fact  into  consideration.  The  random- 
ness arises  out  of  the  hills,  structures,  vegetation,  or  ocean  waves  that  exist. 
When  electromagnetic  waves  are  incident  on  these  surfaces  they  tend  to  diffract 
and  the  resulting  field  is  spread  over  many  scattered  directions. 
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If  we  consider  the  specific  case  of  an  Air  Traffic  Control  (ATC)  sensor 
attempting  to  locate  a beacon  equipped  aircraft,  then  as  a result  of  the  diffuse 
multipath,  the  return  field  viewed  at  the  sensor  resulting  from  a beacon  reply, 
is  a combination  of  a direct  signal,  a specular  multipath,  and  a diffuse  multi- 
path.  The  multipath  signals,  since  they  are  not  coming  from  the  direction  of 
the  aircraft,  may  tend  to  deteriorate  the  performance  of  the  sensor  in  its 
ability  to  determine  the  location  of  the  aircraft. 

The  basic  problem  then  is  for  a fixed  receiving  aperture  and,  with 
some  distribution  of  sensors  over  that  aperture,  how  greatly  does  diffuse 
multipath  degrade  the  angle  estimation  performance?  For  example,  the 
aperture  may  be  a planar  array  of  dipoles  with  a given  horizontal  and  vertical 
distribution.  The  angle  estimation  in  such  an  array  may  be  done  by  some 
form  of  beam  forming  algorithms.  In  contrast  the  array  may  be  a beam 
formed  sum  difference  monopulse  antenna.  For  both  of  these  cases  performance 
evaluations  are  sought. 

To  determine  performance,  it  is  first  necessary  to  characterize  the 
diffuse  multipath  which  in  turn  leads  to  the  characterization  of  the  scattering 
surface.  The  simplest  characterization  is  that  of  a sphere  whose  surface  is 
randomly  perturbed  with  some  statistical  regularity.  The  nature  of  the 
statistical  behavior  of  the  surface  is  of  some  importance  in  determining  the 
structure  of  the  diffuse  field.  The  characterizations  of  the  diffuse  multipath 
can  be  given  in  terms  of  a random  electromagnetic  field  resulting  from  a 
scattering  from  the  random  surface. 

The  problem  discussed  in  this  report  is  an  extension  of  the  work 
reported  elsewhere  to  the  areas  mentioned  above  and  in  so  doing,  it  completes 
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the  study  of  diffuse  multipath  begun  earlier.  To  place  our  approach  in  perspec- 
tive we  will  first  review  the  previous  literature  and  then  briefly  outline  the 
remainder  of  the  report. 

One  of  the  first  analytical  and  experimental  studies  of  scattering  was 
performed  by  Kerr  (1951)  in  an  attempt  to  analyze  the  effect  of  sea  surface 
scattering  on  radar  performance.  The  approach  used  by  Kerr  was  to  obtain 
a scattering  cross  section  for  the  sea  surface.  A more  detailed  approach  was 
undertaken  by  Rice  (1951  ) when  he  modeled  the  surface  by  means  of  a deter- 
ministic polynomial  and  then  solved  the  stochastic  problem.  His  approach  was 
a combination  of  both  deterministic  and  stochastic  analysis  and  directed  at  sea 
surface  scattering.  A different  approach  using  the  Kirchoff  approximation  was 
presented  by  Eckart  (1953)  and  was  directed  at  acoustic  scattering  from  rough 
surfaces.  Equations  for  the  excess  pressure  field  were  obtained  and  then 
averaged.  The  resulting  equations  were  then  solved  for  the  required  quantities. 
This  was  later  followed  by  Ament  (1953,  1956),  using  a similar  approach. 

Using  the  Kirchoff  technique,  Hoffman  (1955)  extended  the  analysis  to 
the  vector  problem  encountered  in  electromagnetic  field  problems.  Other 
attempts  were  made  by  Clark  and  Hendry  (1964)  to  evaluate  backscattering 
effects.  The  prime  reference  using  the  Kirchoff  method  is  Beckmann  and 
Spizzichino  (1963).  Their  results  give  the  field  decomposition  in  terms  of 
different  scattered  directions.  They  are  based  on  a Gaussian  surface  behavior 
with  knowledge  of  surface  standard  deviations  and  correlation  lengths.  These 
approaches  use  the  technique  of  averaging  the  propagation  equations  and  then 
solving.  However,  all  such  techniques  require  some  form  of  closure  argument 
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on  higher  order  moments.  Experiments  yielding  values  for  the  surface  properties 
have  been  discussed  by  Fung  and  Moore,  and  Hayre  and  Moore. 

Most  of  these  previous  results  are  not  suitable  for  an  analysis  of  esti- 
mator performance,  however.  It  was  shown  by  McGarty  (1974),  for  the  case 
of  angle  of  arrival  estimation  with  an  array,  that  a function  called  the  channel 
spread  function  was  necessary  to  evaluate  performance.  This  function  gives 
the  intensity  of  radiation  incident  at  the  receiver  from  all  directions.  This 
function  cannot  be  obtained  directly  from  the  Beckmann  and  Spizzichino  model. 
However,  analyses  by  Kodis  (1966)  and  Barrick  (1968)  provide  exactly  what  is 
necessary  for  the  performance  determination.  Their  analyses  obtain  a scattering 
cross  section  that  relates  the  scattered  radiation  to  a specific  direction  for 
a specified  incident  direction.  The  method  of  analysis  differs  greatly  from 
others  in  that  it  solves  the  scattering  problem  first,  using  a stationary  phase 
technique  (similar  to  Twersky  (1957)),  and  then  an  averaging  technique.  The 
averaging  technique  yields  a very  physical  formula  for  the  scattering  cross 
section.  It  is  in  terms  of  the  average  number  of  scatterers,  the  average 
surface  curvature  (following  Longuet-Higgins ),  and  reflection  coefficients. 

1.  2 Purpose  of  Study 

This  study  has  four  specific  purposes. 

1.  MODEL  DIFFUSE  MULTIPATH  FOR  A SPHERICAL  EARTH. 

To  perform  this  task  it  was  first  necessary  to  characterize  the 
scattering  surface.  Using  both  experimentally  determined  data 
and  analytically  evaluated  models,  the  scattering  surface  can  be 
modeled  as  a Gaussian  random  field.  With  this  model,  it  is 
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then  possible  to  apply  the  Kodis-Barrick  theory  to  the  surface 
and  determine  the  scattering  cross  section.  The  channel  spread 
function  (wave  number  spectrum)  can  be  obtained  for  either  a 
spherical  earth  or  a flat  earth  once  the  scattering  cross  section  is 
evaluated.  The  channel  spread  function  (CSF)  provides  the  power 
density  of  the  diffuse  field  from  every  azimuth  elevation  direction 
at  a given  point. 

2.  EVALUATE  THE  EFFECTS  OF  DIFFUSE  MULTIPATH  ON  THE 
ARRIVAL  ANGLE  ESTIMATION  PERFORMANCE  FOR  DISCRETE 
ELEMENT  ARRAY  ANTENNAS. 

The  characterization  of  the  diffuse  field  provides  the  basis  for 
evaluating  the  degradation  or  improvement  of  angle  estimation 
accuracy.  This  report  considers  only  discrete  element  array 
antennas  (arrays  of  dipoles)  and  evaluates  the  effects  of  diffuse 
multipath  on  the  performance  of  optimum  processors.  By 
optimum  processors  we  mean  the  class  of  maximum  likelihood 
estimators.  Depending  on  the  nature  of  characterization  of  the 
direct  (free  space)  signal,  an  optimum  processor  can  be  obtained. 
These  processors  are,  in  general,  beam  formers  that  search 
over  all  possible  incident  directions  for  the  one  maximizing  the 
processor  output.  The  performance  of  such  processors  (via 
the  Cramer-Rao  [CR]  bound)  acts  as  upper  bounds  to  the  perfor- 
mance of  suboptimum  processors. 
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3.  EVALUATE  THE  EFFECTS  OF  DIFFUSE  MULTIPATH  ON 


MONOPULSE  AZIMUTH  ESTIMATION. 

The  type  of  sensor  chosen  for  the  proposed  DABS  system  is  a 
monopulse  off  boresight  angle  estimator.  Its  realization  may  be 
in  many  forms;  however  one  specific  realization  is  with  a planar 
array  forming  a sum  and  difference  beam.  A great  deal  of 
previous  work  has  been  done  evaluating  the  effects  of  specular 
multipath  on  such  a processor.  However,  only  some  fragmented 
analyses  have  been  performed  on  the  diffuse  multipath  effects. 
Thus  in  this  report  we  have  set  forth  to  evaluate  the  performance 
of  a monopulse  processor  using  the  CSF  developed.  This  proces- 
sor can  then  be  compared  to  the  optimum  processors  to  see  how 
much  of  a margin  of  improvement  can  be  obtained, 

4.  DETERMINE  METHODS  TO  REDUCE  THE  DELETERIOUS 
EFFECTS  OF  DIFFUSE  MULTIPATH. 

In  those  cases  where  diffuse  multipath  causes  degradation  in 
performance,  it  is  important  to  know  the  measures  that  can  be 
taken  to  eliminate  these  effects.  In  particular,  with  monopulse 
azimuth  estimates,  large  vertical  apertures  have  been  suggested 
as  a way  to  eliminate  both  specular  and  diffuse  errors.  It  has 
been  shown  that  for  specular  multipath,  a large  vertical  aperture 
does  eliminate  the  multipath  induced  bias  errors.  Similar  effects 
can  be  studied  with  monopulse  and  diffuse  multipath  using  the  CSF 
and  the  performance  evaluations  obtained. 
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These  are  the  main  purposes  of  this  study.  Of  particular  interest 
is  the  geometry  that  one  finds  in  an  air  to  ground  surveillance  link  with 
special  emphasis  on  the  DABS  implementation.  The  results  obtained  are 
more  general,  however,  and  thus  this  study  addresses  many  other  issues 
not  specifically  outlined  in  this  section. 

1.  3 Results  of  the  Study 

This  report  has  provided  many  results  that  apply  directly  to  DABS  and 
the  general  ATC  problem  as  well  as  to  areas  of  a more  general  nature.  The 
following  are  the  specific  results  and  conclusions  reached. 

1.  DIFFUSE  MULTIPATH  CAN  BE  COMPLETELY  CHARACTERIZED 
BY  THE  CHANNEL  SPREAD  FUNCTION  FOR  A WIDE  CLASS  OF 
SURFACE  CHARACTERIZATIONS. 

As  previously  discussed,  the  diffuse  field  is  dependent  upon  the 
surface  characterization.  In  this  report  a rather  robust  set  of 
surface  characterizations  has  been  chosen,  and  a statistical 
description  of  the  scattered  electromagnetic  field  has  been  obtained. 
The  second  order  statistics  of  the  field  are  obtained  by  considering 
the  field  as  received  at  the  aperture  of  a sensor  resulting  from 
an  excitation  by  a source  (a  beacon  equipped  aircraft)  at  a loca- 
tion well  defined  with  respect  to  the  sensor.  The  result  is  the 
channel  spread  function,  which  provides  the  energy  density  at  the 
receiving  aperture  in  terms  of  azimuth  and  elevation.  The  CSF 
depends  upon  the  sensor  and  source  geometries,  the  surface 
statistics  in  terms  of  its  r.  m.  s slope,  the  source  antenna  pattern. 
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as  well  as  shadowing  effects.  The  CSF  provides  a complete 
analytical  characterization  of  the  diffuse  field. 

2.  DIFFUSE  MULTIPATH  CAN  BE  SIGNIFICANTLY  REDUCED  BY 
SHADOWING. 

Shadowing  results  when  surfaces  are  not  directly  irradiated  by 
the  incident  electromagnetic  field  due  to  blockage.  The  charac- 
terization of  this  type  of  effect  is  quite  difficult  for  random  sur- 
faces, but  there  exist  several  viable  models,  particularly  the  one 
by  Wagner.  Using  this  model,  the  total  power  of  the  diffuse 
multipath  is  decreased  significantly  at  low  elevation  angles  as  a 
result  of  surface  blockage.  The  difference  between  including 
shadowing  effects  and  disregarding  them  at  angles  on  the  order 
of  1 to  2^  is  on  the  order  of  1 0 to  20  dB  in  power.  Thus  when 
shadowing  is  included  in  the  model  as  it  is  in  the  actual  propagation 
medium,  the  total  diffuse  power  can  be  significantly  reduced. 

OPTIMUM  PROCESSORS  UTILIZE  THE  DIFFUSE  COMPONENT 
TO  IMPROVE  ANGLE  ESTIMATION  PERFORMANCE. 

The  optimum  processor  is  a maximum  likelihood  estimator  of 
the  arrival  angle.  To  obtain  an  estimate  with  such  a processor, 
it  is  necessary  to  form  many  beams  and  determine  on  which  one 
the  signal  is  arriving.  This  entails  a great  deal  of  processing  and 
uses  all  the  information  presented  the  sensor  over  its  entire 
aperture.  Such  processors  require  outputs  from  each  dipole  and 
as  such  are  quite  complicated.  However  when  such  processors 
are  used  in  the  presence  of  diffuse  multipath,  which  is  centered 
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at  the  azimuth  direction  of  the  incident  plane  wave,  these  proces- 
sors can  utilize  this  information  to  improve  the  estimate.  In 
effect  the  optimum  processor  utilizes  the  information  regarding 
source  location  in  the  diffuse  field  that  generated  it,  and  reduces 
source  location  uncertainty.  The  result  is  presented  in  terms 
of  the  Cramer-Rao  bound  on  angle  estimate  variance.  For  the 
optimum  processor,  a decrease  in  the  estimate  variance  can, 
in  theory,  be  obtained  by  processing  both  the  direct  and  diffuse 
fields. 

4.  MONOPULSE  AZIMUTH  PERFORMANCE  IS  ALWAYS  DEGRADED 
BY  THE  PRESENCE  OF  DIFFUSE  MULTIPATH. 

In  contrast  to  an  optimum  multibeam  processor,  the  two  beam 
sum-difference  monopulse  processors  always  have  their  perfor- 
mance degraded  by  diffuse  multipath.  This  can  best  be  understood 
by  noting  that  for  a monopulse  processor,  both  beams  point 
toward  the  target  whose  position  is  to  be  estimated.  The  result 
is  that  the  diffuse  component  acts  as  noise  and  raises  the  effec- 
tive noise  level.  The  azimuth  estimation  performance  is  then 
degraded.  The  degradation  can  be  a factor  of  2 to  3 above  the 
noise  standard  deviation  and  almost  an  order  of  magnitude  above 
the  performance  obtained  with  an  optimum  processor.  However, 
for  the  purposes  of  a DABS/ATC  sensor,  the  performance  is 
still  acceptable  even  in  the  worst  diffuse  multipath  conditions. 
Moreover,  the  diffuse  errors  are  still  an  order  of  magnitude 
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smaller  than  those  generated  by  specular  multipath.  An  important 
point  to  note  however  is  that  in  contrast  to  noise,  the  diffuse 
multipath  azimuth  error  does  not  depend  on  the  signal  strength. 

This  is  because  as  this  signal  level  is  raised,  the  diffuse  level 
is  also  raised  proportionately. 

5.  VERTICAL  APERTURE  IS  USEFUL  IN  REDUCING  MONOPULSE 
DIFFUSE  MULTIPATH  AZIMUTH  ERRORS. 

The  CSF  shows  that  the  diffuse  multipath  is  concentrated  at  the 
horizon  corresponding  to  small  negative  elevation  angles.  This 
implies  that  for  targets  at  angles  above  1^  to  2^,  vertical  aperture 
can  be  used  to  reduce  the  effect  of  diffuse  multipath.  At  smaller 
elevation  angles,  the  diffuse  field  is  self  limiting  as  a result  of 
shadowing  and  spreading  in  azimuth.  The  net  result  is  that  for 
angles  greater  than  1.  5°  vertical  aperture  helps,  and  for  angles 
less  than  this  diffuse  multipath  is  itself  small.  Thus  vertical 
aperture  works  best  where  the  diffuse  field  is  strongest.  This 
fact  is  demonstrated  with  data  from  radial  flights  performed  at 
DABS  Experimental  Facility  (DABSEF). 

These  five  results  are  the  main  conclusions  of  this  study.  They 
show  that  the  diffuse  field  can  be  characterized  and  its  effect  on  azimuth 
accuracy  analytically  evaluated. 

1.  4 Summary  of  the  Report 

The  main  text  of  this  report  presents  a detailed  summary  of  the  analysis 
of  the  diffuse  field  and  the  array  performance.  The  specifics  of  the  analyses 
are  contained  in  the  Appendices. 
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In  Section  2,  the  model  for  the  received  signal  is  developed.  It  is 
assumed  that  future  analyses  will  be  concerned  with  arrays  so  that  the  spatial 
model  is  for  point  detectors  over  a given  aperture.  The  time  behavior  is  also 
defined  in  terms  of  the  bandpass  nature  of  the  received  signal.  The  model  of 
the  diffuse  signal  is  presented  and  its  statistics  discussed.  The  resulting 
voltage  output  is  derived  and  related  to  the  incident  electric  field.  The  main 
result  of  this  modeling  is  the  introduction  of  the  channel  spread  function. 

Section  3 will  evaluate  the  channel  spread  function  for  a spherical  earth. 
The  surface  irregularities  are  modeled  as  a homogeneous  time  invariant 
Gaussian  random  field.  Using  the  Kodis -Barrick  theory,  the  scattering  cross 
section  is  evaluated.  Shadowing  effects  are  included  in  the  model  using  the 
results  of  Wagner,  Specific  results  on  the  evaluation  of  the  CSF  for  various 
geometries  and  surface  roughness  are  presented. 

The  performance  and  structure  of  various  angle  estimators  are 
discussed  in  Section  4,  The  signal  model  developed  in  Section  2 is  used  as 
a basis  for  various  estimator  structures.  The  optimum  maximum  likelihood 
estimators  for  both  spread  and  unspread  channels  are  discussed  in  this  section. 
A detailed  comparison  of  structure  and  performance  is  made  under  the  varying 
assumptions  of  signal  structure.  Particular  emphasis  is  placed  upon  the  mono- 
pulse processor,  and  its  performance  is  compared  to  that  of  the  optimum 
processor. 
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SECTION  2 


THE  SIGNAL  MODEL 


When  a signal  is  transmitted  from  a source  at  a given  location  and 
received  by  a sensor  at  another  pointy  the  total  signal  can  be  broken  down  into 
four  parts.  They  are  the  free  space  portion,  the  specular  multipath  portion, 
the  diffuse  multipath  term,  and  noise.  The  specific  description  of  the  received 
signal  depends  upon  the  source  characteristics,  the  propagation  nature  of  the 
medium,  and  the  type  of  sensor  used  to  measure  the  fields  incident  upon  it. 

To  allow  the  extension  of  the  results  obtained  in  this  report  to  as  wide  a 
variety  of  sensors  as  possible,  the  received  signal  is  best  described  as  re- 
ceived by  a point  detector  (a  single  dipole)  located  at  a specific  point  in  the 
receiving  aperture.  With  this  type  of  description,  an  extension  to  arbitrary 
sensor  beam  patterns  can  be  made  by  appropriate  phasing  and  weighting  of  the 
individual  point  sensors. 

The  aperture  of  the  receiving  sensor  can  be  defined  with  respect  to  an 
arbitrary  coordinate  system  as  shown  in  Fig.  1,  In  this  aperture,  Ji  , we  can 
locate  a specific  point  defined  by  a vector  r.  If  a sensor  is  placed  at  this 
point,  then  its  output  voltage  as  a function  of  time  is  given  by  r(t).  This 
voltage  depends  upon  the  incident  electromagnetic  field  that  is  propagating  at 
a frequency,  f^,  the  carrier  frequency.  Thus  r(t)  is  actually  a signal  modulated 
at  this  carrier  frequency.  In  general  the  modulation  on  r(t)  is  narrowband 
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Fig.  1,  Geometry  of  receiving  aperture  . 
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with  respect  to  so  that  r(t)  can  be  represented  in  what  is  termed  as  complex 
envelope  form.  Namely 

r(t)  = y/z  [r(t)  exp  (jZirf^t)]  (2,  1 ) 

where 

r(t)  is  the  complex  envelope  or  r(t). 
r(t)  is  a slowly  varying  time  function  (e.  g.  , a pulse)  that  represents 
the  baseband  modulation  of  the  carrier.  It  is  therefore  sufficient  to  deal  with 
only  the  complex  envelope,  r(t). 

The  first  portion  of  the  received  signal,  the  free  space  portion,  is 
given  in  the  following. 

Sr(t)  = f(t  - R/c)b  exp  • £)  (2.  2) 

is  the  energy  in  the  signal  when  transmitted 
is  the  complex  envelope  of  the  time  variation 
is  the  complex  amplitude  and  phase  variations  due  to  path 
propagation 

is  the  range  from  the  source  to  the  sensor 
is  the  propagation  velocity 
is  the  wave  number 

is  the  direction  cosine  vector  indicating  the  incident  direction 
of  the  wave  front 

The  time  variation  has  the  property  that  over  an  interval  (0,  T),  it  is 
normalized  such  that 


where 


^0 

f(t) 

fsj 

b 


R 

c 
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T 


f(t)  f * (t)dt  = 1 


(2.  3) 


where 


represents  complex  conjugate 


The  second  portion  of  the  received  signal  is  the  specular  multipath 
return  and  is  given  by- 


having  been  reflected  from  the  earth's  surface  and  also  possibly 
being  altered  by  an  antenna  gain  change  at  the  source 
R*  is  the  path  length  over  the  new  path 

jj*  is  the  direction  cosine  vector  of  the  multipath  wave. 

The  diffuse  signal  can  also  be  characterized  using  this  format. 


However,  it  must  first  be  noted  that  the  diffuse  field  is  assumed  to  come  from 
a continuum  of  directions  and  not  just  one  as  the  free  space  and  specular 
returns.  To  represent  this  continuum,  we  use  an  integral  representation  in 
the  following  form 


(2.4) 


where 


p is  the  complex  reflection  coefficient  that  results  from  the  wave 
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where 


R(p')  now  represents  the  path  length  from  source  to  sensor  over  the 

diffuse  multipath  source  arriving. at  the  sensor  from  direction  p' 
b(M*)  is  the  amplitude  and  phase  dependence  of  that  particular  direction. 

To  fully  characterize  the  diffuse  term,  a complete  description  of 
b (p)  must  be  made.  Since  the  reflecting  surface  giving  rise  to  the  surface  is 
random,  then  so  too  is  b (p).  Furthermore,  since  the  field  coming  from 
direction  p may  most  likely  arise  from  many  points  independently  radiating  in 
that  direction  (as  will  be  shown  in  the  next  section),  then  using  standard 
limiting  arguments,  one  can  show  that  b (p)  is  a complex  Gaussian  random 
variable  for  each  p.  Furthermore,  it  can  also  be  shown  to  be  zero  mean. 

To  fully  characterize  it,  all  that  is  needed  is  the  covariance.  Experience 
indicates  that  spatially  the  diffuse  field  is  a homogeneous  random  field  (see 
Capon)  so  that  for  the  spectral  (wave  number)  portion,  we  obtain  the  following. 

For  a source  whose  free  space  direction  is  p,  we  have 

E[b  « (m")]=  K(p,  - m")  (2.  6) 


where 

6(p)  is  the  delta  function. 

This  covariance  states  that  radiation  coming  from  different  directions  is 
independent.  The  term  K(p>  p*  )represents  the  spectral  density  from  direction 
p*  for  a source  at  p.  This  function  is  called  the  channel  spread  function 
(CSF)  because  it  represents  the  spreading  of  the  signal  due  to  the  diffuse  field. 


16 


The  fourth  and  final  term  in  the  received  signal  is  the  noise,  a zero 
mean  narrow  band  white  Gaussian  noise  process  w(t).  For  this  process  we 
have 

E[w(th^^  * (t  + t)]  = Nq  6(t)  (2.7) 

Thus  the  received  signal  at  point  r is  given  by 

?(t)  = ^^(t)  + Sg(t)  + Sj^(t)  + w(t)  (2.  8) 


It  is  often  more  convenient  to  deal  with  the  demodulated  form  of  r(t) 
at  each  sensor.  This  is  obtained  by 


•T+R/c 


r = 


r(t)f  (t  - R/c)dt 


'R/c 


(2.9) 


where 

r(t)  is  demodulated  by  the  delayed  complex  envelope  of  the  signal. 

The  result  will  be  to  eliminate  the  envelope  variations  in  both  the  direct  and 
specular  paths  and  introduce  a known  term  in  the  diffuse  path.  The 
demodulated  diffuse  term  is 


where 


g(/j')  b(/j')  exp  • r)dfx' 


(2.  10) 


T+R/c 


g(M')  = 


f(t  - R(iU')/c)r'''  (t  - R/c)dt 


(2.  11) 


Now  over  most  viewing  directions,  g(|u’)is  almost  unity  by  the  nature  of  f(t) 
so  it  can  be  neglected. 

We  can  now  relate  this  signal  model  to  the  field  incident  on  the  sensor. 

Let  us  first  consider  the  diffuse  field  only.  Let  f(t)  dE(0,4>)  be  the  field 

coming  from  the  solid  angle  0,  0 + d0,  <|),  <|)  + d<^  at  time  t.  Since  this  represents 

the  diffuse  field  it  is  . random.  Assume  that  the  sensor  is  a dipole  with  cross 
2 

section  X /4tt,  where  X is  the  wavelength.  Then  in  terms  of  the  electric  field, 
the  mean  square  voltage  output  can  be  written  as 


= Eq  f(t)f  (t) 


dE(0,cj,)|^Jf(t)f  * (t)T 


(2.  12) 


where  we  have  equated  it  to  the  mean  square  of  the  voltage  output  using  the 
diffuse  model  developed  in  this  section,  and  define  0 and  (f>  as  shown  in  Fig,  2, 
Now  if  we  also  assume  the  electric  field  to  be  independent  from  direction  to 
direction,  we  define 


E[dE(0,4,)  dEt  (0',4,')1  = 


K(0,  4))d0d4>; 


0/0',  4./4,' 
0 = 0',  4>  = 4>' 


(2.  13) 


where 

t is  the  complex  conjugate  transponse, 

^ 2 2 

The  units  of  K(0,  (f))  are  in  volts  /m  -steradian.  Using  this  definition  and 
equating  terms  in  the  mean  square  voltage  output  we  find  that 
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Fig.  2.  Diffuse  multipath  geometry. 
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(2.  14) 


K(M.iu') 


^2^ 
X T 

^firET 


K(Q,<^) 


Thus  if  from  the  geometry  and  propagation  equations  we  can  find  the  electric 
field  and  its  statistics,  then  K(0,  cj))  is  obtained  and  concurrently  K(p,  jj*).  It 
is  the  purpose  of  the  next  section  to  show  how  this  term  is  obtained. 

In  general,  when  dealing  with  arrays,  the  total  demodulated  received 
signal  is  represented  by  the  output  of  N sensors  located  at  positions  r^, 
i = 1,  . . . , N,  If  we  introduce  an  N x 1 vector  fii(p),  called  the  delay  vector, 
given  by 


exp  (jk^M  • r ^ ) 


_exp  (jkpM  • rj^)_ 


(2.  15) 


then  we  can  easily  write  the  total  output  of  an  array  in  vector  form.  Specifi- 
cally, if  we  let  r^  be  the  demodulated  output  at  location  r.,  then  define  r as  the 
N X 1 vector  of  such  outputs  then 
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where 


w is  now  an  N x 1 zero  mean  Gaussian  random  vector  with  covariance 
Kw 

With  this  model  for  a received  signal,  one  can  now  ask  how  to  estimate 
the  direction  cosines /it  , or  equivalently  the  azimuths  and  elevations.  For  the 
purposes  of  this  report  we  are  interested  in  only  azimuth  so  that  fj,  will  be  a 
scalar.  However,  the  general  or  canonical  form  of  the  received  signal  is 
still  of  use. 

The  model  we  have  developed  has  provided  the  following. 

1.  A description  of  the  received  signal  in  terms  of  a free  space  path, 
a specular  multipath,  a diffuse  multipath,  and  noise 

2.  A time  demodulated  waveform 

3.  A model  for  the  diffuse  component  that  shows  it  to  be  a zero  mean 
homogeneous  Gaussian  random  field 

4.  A relationship  between  the  model  and  the  electromagnetic  field 
incident  on  a dipole  yielding  the  measured  output  voltage 

5.  A canonical  form  for  an  array  output  useful  for  estimating 
arrival  angles. 

With  this  model,  in  the  next  section  we  shall  analytically  obtain  K(iu,)u’)  from 
first  principles,  and  in  Section  4 evaluate  various  processors  to  estimate  0 
the  azimuth,  or  fj,  the  direction  cosine. 
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SECTION  3 


THE  CHANNEL  SPREAD  FUNCTION 

The  channel  spread  function  is  a measure  of  how  the  diffusely  scattered 
radiation  is  distributed  in  azimuth  and  elevation  with  respect  to  the  aperture 
of  a receiving  antenna.  In  this  section  we  shall  derive  this  expression  in  some 
detail  showing  how  it  arises  from  the  random  nature  of  the  surface  and  the 
laws  of  propagation.  In  order  to  solve  the  equations,  certain  assumptions  must 
be  made  and  the  implications  and  limitations  of  these  assumptions  will  be 
discussed. 

In  this  section  we  first  obtain  the  CSF  in  terms  of  a certain  set  of 
parameters  and  show  how  it  arises  naturally  as  a result  of  surface  scattering. 
The  parameters  needed  in  the  definition  are  the  scattering  cross  section  per 
unit  area,  and  suitable  Jacobian  and  shadowing  effects.  In  the  remainder  of 
this  section  these  effects  are  evaluated. 

3.  1 The  CSF 

In  the  previous  section  we  evaluated  the  relationship  between  the  model 
of  the  received  voltage  at  the  output  of  a dipole  and  the  electric  field  incident 
on  the  dipole.  This  field  was  described  as  having  come  from  a continuum  of 
azimuth  (9)  and  elevation  (cj>)  directions  with  respect  to  the  dipole.  Consider 
the  geometry  in  Fig.  3.  A source  of  radiation  is  at  point  T and  it  is 
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dA 


Fig.  3.  Geometry  of  incremental  scattering  surface  . 
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received  at  point  S.  The  source  has  power  P and  an  antenna  with  gain 

s 

G^(0^,  <j)^)  where  0^  and  (j)^  are  the  azimuth  and  elevation,  respectively,  with 
respect  to  the  radiated  energy.  These  two  points  are  located  over  a spherical 
earth  of  radius  R„,  The  transmitter  is  at  a height  z , and  the  sensor  at 
height  Now  consider  radiation  coming  to  the  sensor  S from  solid  angle 

(j),  (j)  + d<j)  and  0,  0 + d0.  In  the  figure  we  show  that  if  we  let  $"=(})+  tt/2,  then 
this  solid  angle  generates  a surface  area  dA  from  which  power  must  be  scattered 
so  as  to  be  received  from  this  direction.  Let  s^  be  a vector  from  the  sensor 
at  S to  the  center  of  dA,  and  s^  a vector  from  T to  the  center  of  dA.  A 
coordinate  system  is  centered  below  S on  the  surface  of  the  earth.  Also  assume 
the  radiation  has  wavelength  X, 

Let  0.  and  (j>.  represent  the  azimuth  and  elevation  at  the  surface  dA 

2 

to  the  transmitter  at  T.  The  transmitter  radiates  a power  density  (watts/m  ) 

2 

of  (Pq/4tt  R ) G,p(0^,  <j)^)  at  a range  R and  in  direction  0^,  cj>^.  Then  the  power 
density  on  dA  is  (Pq/^ttI^^I  G^(0^,  4>^).  The  surface  has  area  dA  so  that  the 
total  power  can  be  evaluated.  Let  represent  the  wave  vector  defining  the 
direction  of  radiation  coming  from  T at  dA,  and  let  k represent  the  wave 

— ^5 

vector  defining  propagation  from  dA  to  S.  To  find  how  much  of  the  incident 
power  is  scattered  from  k^  to  on  dA,  namely  from  transmitter  to  receiver, 
we  need  to  define  a scattering  cross  section,  cr  (k  , k. ).  This  scattering  cross 
section  reflects  how  effective  dA  is  in  transmitting  power  from  the  incident 
direction  into  the  scattered  direction. 

Having  the  power  scattered  toward  the  receiver  S,  it  is  then  possible 
to  obtain  the  power  density  (watts/m  ) at  S as  follows 
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Pf)  G (0.,  (j). ) 2 

2"j O' (k^,  ^ ) dA  (watts/m  ) 


(3.  1) 


To  obtain  K(0,  cj)),  the  mean  square  voltage  density  (volts  /m  -steradian), 
we  use  the  radiation  resistance  of  120tt  ohms  (see  Friis  and  Lewis)  and  the 
solid  angle,  dOdcj).  Thus 


Pq 

♦i)  “-ti-hs)  dA 

(4Trr|sJ 

r\ 

1 1 2 1 ZOtt  d0d(j) 

1 s^l 

(3.2) 


where  0 and  (j>  are  the  azimuth  and  elevation  of  dA  with  respect  to  S. 

Once  K(0,  (j))  is  obtained  the  CSF,  K(^),  follows  directly  from  the 
evaluation  in  Section  2.  It  can  be  noted  that  in  the  definition  of  K(0,  cj)),  the 
first  term  depends  on  the  geometry  of  the  two  sensors,  the  second  term  on 
the  nature  of  the  scattering  surface,  and  the  third  term,  a Jacobian,  on 
the  relationship  of  the  scattering  surface  mean  value  to  the  viewing  angles. 
Each  of  these  terms  will  now  be  evaluated  starting  with  the  scattering  cross 
section. 

3^  2 Scattering  Cross  Section 

Consider  the  geometry  as  described  above  and  let  ^ vector 

pointing  from  S to  dA.  If  we  define  an  arbitrary  coordinate  system  centered 
in  dA,  then  by  using  the  Stratton- Chu  formulation,  we  can  obtain  the  field 
from  this  surface  area  at  S.  Let  this  field  be  given  by  E(r^)  where  r^  is  a 
vector  from  the  center  of  dA  to  S that  is  in  the  same  direction  as  k . Now 
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assume  that  the  field  incident  on  dA  is  a plane  wave  with  a wave  vector  k. 

with  respect  to  the  coordinate  system  defined  at  dA.  Let  the  amplitude 

of  the  field  be  and  assume  it  has  an  arbitrary. polarization.  Then  the 

2 

scattering  cross  section  (in  m ) is  defined  as 


0-  (k;,  k ) = 


4ttR 


(3.  3) 


where  R is  the  distance  from  S to  dA,  and  E [ ] represents  the  expected 
value  (see  Kerr,  p.  33). 

The  expectation  is  necessary  because  the  surface  is  irregular  and 
random  in  nature,  so  that  the  field  scattered  from  it  is  random.  Let  ^(x,  y) 
be  the  height  of  the  surface  at  a specific  x,  y coordinate.  For  this  analysis 
we  assume  that  ^(x,  y)  is  a Guassian  random  variable  with  zero  mean 

E[e(x,y)]=0  (3.4) 

and  a covariance  of  the  form 


E[^(x,  y)  ^(x',  y')]  = exp  (-r^/cr^) 

where 

cTj^  is  the  r.  m.  s.  surface  height 
0*^  is  the  correlation  length 

also 

7 2 * ? 

r^  = (X  - X* ) + (y  - y r 


(3.5) 


(3.  6) 


26 


Thus  over  any  incrementally  small  area  dA,  this  planar  description  of  the 
surface  height  holds  for  a coordinate  system  centered  in  dA,  The  correlation 
function  for  the  surface  has  been  chosen  to  satisfy  certain  regularity  conditions 
mentioned  by  Bar  rick  [Z], 

In  order  to  solve  the  boundary  value  problem  that  arises  from  the 
Stratton-Chu  formulations,  we  assume  that  the  surface  is  a perfect  conductor. 
Then,  as  shown  in  Appendix  I,  the  resulting  integral  over  the  surface  area  is 
in  a form  that  can  be  solved  if  stationary  point  approximations  are  valid. 

These  approximations  mean  that  the  average  radii  of  curvature  of  the  surface 
are  much  larger  than  a wavelength  so  that  points  that  lead  to  scattering  do 
so  in  a specular  manner.  Thus  with  this  Approximation,  we  have  a specular 
point  theory  valid  for  high  frequencies,  generally  at  UHF  and  above. 

The  net  result  of  this  approximation  is  that  (k^,  k^)  can  be  written  as 

) = un  R^(k.,k  ) (3.7) 

—X  ““S  X ^”X  ^“S 

where 

n^  = the  average  number  of  specular  reflectors 

= the  average  of  the  product  of  the  radii  of  curvature 

2 

R (k.,  = the  reflection  coefficient  giving  the  depolarization  loss 

from  a tilted  surface. 

If  we  let  dA  become  incrementally  small,  then  we  can  equivalently  define  the 
scattering  cross  section  per  unit  area  as 
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(3.8) 


0-  (k., 


!ls>  = 


where 

= the  average  number  of  specular  points  per  unit  area 

The  average  of  both  the  specular  points  and  the  product  of  the  radii 
of  curvature  can  be  obtained  using  the  hypothesized  statistical  distribution  of 
the  surface.  This  has  been  done  by  Barrick  [1]  and  is 
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where 

2 

s 


and 


the  square  of  the  rms  surface  slope 


tan  7 


I 2 ! 2 

\ sin  (j).  - 2 sin  sin  cos  0^  + sin 

cos  (|).  + cos  c|) 


(3.  10) 


where 

7 is  the  angle  between  the  local  normal  on  dA  and  the  local  surface 
normal  of  the  specular  point. 

The  angles  and  are  the  angles  between  the  local  normal  on  dA  and  the 

propagation  vectors  k.  and  k . 0 is  the  angle  of  scattering  in  the  local  x,  y 

*“~S  s 

plane  between  k.  and  k (see  Fig.  4). 

— 1 — s 
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Fig.  4.  Spherical  earth  scattering  geometry. 
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Note  that  n . represents  a density  of  specular  scatterers  per  unit  area 
and  it  depends  upon  the  location  of  the  surface  and  the  statistical  characteristics 
of  the  surface.  For  surfaces  with  large  rms  slopes,  s,  n^  can  tend  to  be 
large. 

Barrick  has  also  evaluated  the  average  product  of  the  radii  of  curvature. 

It  is 


2 4 

^ = 0.  1378tt  (T£  sec  7 

"~2 

s 


Combining  these  we  obtain  the  scattering  cross  section 


(3.  11) 


R^(k.,k  ) 
—1  —s' 


(3.  12) 


2 

The  result  shows  that  cr  (k*,  kg)  depends  upon  y,  s and  R (k^,  k^),  and 

is  independent  of  (r  ^ and  cr^  separately,  i.  e.  , it  depends  upon  the  rms  surface 

slope  and  it  is  independent  of  wavelength. 

The  wavelength  independence  results  from  the  fact  that  we  used  the 

specular  point  approximation.  Thus  it  is  truly  wavelength  independent  only  if 

2 

Tfis  greater  than  X . From  the  dependence  of  oncr^,  s and  y,  we  see  that 
this  assumption  depends  upon  the  surface  statistics  and  the  geometry  of  fhe 
scattering  point. 
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The  reflection  coefficient  has  been  evaluated  by  Mitzner,  For  horizontal 


to  horizontal  or  vertical  to  vertical  polarization  we  have 


R (k.,  k ) = 

l’  —8 


’ -sin  (j).  sin  ()>  sin  0 + a a 

X S S M j 

— — ^ 

4 sin  § cos  § 


(3.13) 
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^1 
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s 

(3.  14a) 
(3.  14b) 
(3.  14c) 


Similar  results  apply  for  cross  and  circular  polarization. 

The  fundamental  constraint  or  assumption  using  this  method  is  that 
the  reflectors  are  quasi- specular  which  assumes  that  the  wavelength  of  the 
radiation  is  small  with  respect  to  the  radii  of  curvature.  Thus,  the  results 
are  valid  on  the  short  wavelength  limit.  Beckmann  and  Spizzichino  (p.  89  ) 
also  evaluate  a mean  power  reflection  coefficient,  which  when  suitably  nor- 
malized, is  identical  to  (3.  7)  but  without  the  reflection  coefficient  term. 
They  state  that  their  result  is  valid  for  very  rough  surfaces  where 


g = 


2it 


-J-  (cos  ()).  + cos  ())^ 


2 

» 1 . 


(3.  15) 


This  agrees  with  the  limitation  on  the  Rayleigh  criterion  for  rough  surfaces, 
which  states  that  a surface  can  be  considered  smooth  if 
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> 


(3.16) 


^ h 8 cos  (j) 

Thus,  for  low  angle  scattering  (cj>  tt/Z)  we  have  the  constraint  that  the  model 
we  are  considering  is  valid  for 


"7" 


» 


1 


^ = tt/Z  - cj)  . 


(3.  17) 


For  example,  at  3°,  this  means  that  must  be  greater  than  Z.  5. 

Thus,  for  X = 1 foot,  for  example,  the  surface  heighLs  for  diffuse  multipath 
are  on  the  order  of  8 ft  or  more.  Such  heights  are  common  in  hilly  terrain. 
It  furthermore  states  that  at  such  angles  small  rocks,  etc.  , appear  smooth. 

Thus,  the  model  that  we  are  considering  is  valid  under  the  following 
constraints : 

1.  The  radius  of  curvature  is  large  compared  to  a wavelength. 

This  means  that  the  specular  point  approximation  can  be  used. 
This  implies  that  diffraction  can  be  neglected  from  each 
independent  reflector. 

Z.  Multiple  scattering  can  be  neglected.  The  problems 

associated  with  multiple  scattering  are  significant  and  are 
discussed  in  McGarty  [3].  It  basically  relates  to  the 
closure  property  associated  with  solving  for  the  moments  of 
random  integral  equations. 

3.  That  cj)^  is  such  that  the  surface  cannot  be  considered  smooth. 
This  implies  (3.  15)  is  valid.  This  allows  for  the  scatterers 
to  be  added  incoherently. 


3Z 


The  first  assumption  is  the  short  wavelength  approximation,  while 

that  for  (3)  implies  a very  rough  surface.  Now  as  (3.  12)  demonstrates, 

2 

O'  (kji  j5g)  depends  only  on  the  ratio  s and  not  on  or  independently. 

However,  this  dependence  is  implicit  under  assumption  (3), 

We  can  attempt  to  extend  the  model  by  augmenting  the  saddle  point 
integration  with  other  techniques  (see  Felsen  and  Marcuvitz,  Chpt.  4). 

However,  the  intuitive  expres sion  for  the  scattering  cross  section  breaks 
down.  In  conclusion  then,  the  expression  for  the  scattering  cross  section 
should  be  adequate  for  a wide  range  of  surfaces  found  in  operational  environ- 
ments. 

3.  3 Surface  Jacobian 

The  third  portion  of  the  channel  spread  function  relates  the  incremental 
surface  area  to  the  solid  angle  that  generates  it.  This  is  in  effect  a Jacobian 
transforming  surface  area  to  solid  angle.  This  function  can  be  evaluated  for 
various  surfaces  such  as  planar,  spherical  and  others.  The  first  portion  of 
the  CSF  also  depends  upon  the  geometry  with  respect  to  a given  0,c})  viewing 
angle.  In  this  section  we  will  relate  0 and  to  the  sensor  and  source  locations 
and  the  scattering  angles  at  the  specular  points.  Specifically  we  shall  consider 
the  geometry  for  a spherical  earth,  that  for  a flat  earth  having  been  done 
elsewhere  (see  McGarty  [2]). 

Consider  a receiver  as  shown  in  Fig.  4 with  azimuth  0 and  elevation  c|) 

to  a scattering  surface.  Let  n^  be  the  local  normal  on  the  scattering  surface, 

and  k.,  k be  defined  as  before.  Let  the  source  be  at  an  angle  from  the 
— 1 — s 

receiver  referenced  to  the  center  of  the  earth.  The  ground  range  to  the 
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transmitter  is  R = or*  R_.  Let  the  receiver  be  at  height  z and  the  transmitter 

hj  1 

at  z^.  In  this  case  also  assume  that  the  transmitter  and  receiver  are  coplanar. 

Now  let  a be  the  angle  to  the  center  of  the  scattering  surface.  It  is 
readily  shown  that  given  a can  be  obtained  from  the  relationship 

j 

z . + R„(  1 - cosor) 

tan  cb  = 5 1 . (3.  18) 

^ R^  sinor  ' ' 

Thus,  to  obtain  K(0,c()),  we  have  for  each  0,  a defined  incremental  surface 
at  an  angle  or  relative  to  a spherical  coordinate  system  centered  at  the  earth^s 
core.  Similarly,  the  surface  is  also  defined  by  a p given  by 

0 = 2 - P • (3.  19) 

From  the  results  of  the  previous  subsection  we  saw  that  and  0 must 

X s s 

be  defined  for  each  0,c|)  surface.  In  the  Appendix,  we  show  that 

4)s  = 5 - 4)  + a (3.  20) 

and  (j).  and  0^  are  rather  complicated  but  readily  evaluated  functions  of  0 and  (j). 
Similarly  |sj|  and  |s2l  are  given  by 

1^1 1 “ ^ ^^E^^  - cosor)  + z j)  (3.  21) 

“ cos^p  sin^or  + ((z2  + H^^;)  sinor'  - sinp  sinor  R^-)^ 

+ ((Z2  + Rj,)  cosor'  - cosor  R^-)^  • (3.  22) 
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Finally,  from  the  geometry  we  note  that 


2 

dA  = Rj,  sina  da  dp 


(3.  23) 


Using  the  functional  relationships,  dA  can  be  shown  to  reduce  to  (see  Appendix  II) 


dA  = 


111 


^sc^a  - csca  coto^^  + ijj 


dQ  d({) 
sina 


(3.  24) 


Thus  the  spatial  spectrum  of  the  field  is 


"*  ^0  ^®i’  ‘I’i^ 

K(e,4))  = — V,  ■■  ,2  4>i) 


(4tt) 


\lz\ 


sina 


1 - cosall  + 


(3.  25) 


where  cj).,  <\>  , 0 and  a all  are  defined  in  terms  of  0 and 
X s s 

Before  continuing  we  must  consider  the  limiting  factors  on  a and  p. 
Clearly,  restrictions  on  p are  the  same  as  those  on  0;  consequently  they  are 
not  important.  However,  a can  only  be  so  large  since  we  must  consider  the 
effect  of  the  horizon.  Consider  the  geometry  in  Fig.  5.  Here  we  show  the 
largest  value  of  a,  a.,  and  the  corresponding  d>  value.  This  occurs  when  the 
ray  from  the  source  just  touches  the  surface  of  the  earth.  Thus 


R, 


cos  Q.  = -r 

* R_  + z, 


(3.  26) 
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Fig.  5.  Horizon  cutoff. 


36 


is  the  maximum  Thus  all  acceptable  a are  those  such  that 


cos  a > 


R, 


(3.  27) 


which  implies 


1 - 


(3.  28) 


For  the  case  where  z.  <<R„ 

1 E 


or 


♦ 


(3.  29) 


To  find  the  ground  range  to  the  target  at  this  maximum  angle  we  note  that 


cos  5a  = 


R, 


E + ^2 


(3.  30) 


and  the  maximum  ground  range  is 


R 

max 


or* 


(3.  31) 


In  general,  z^<<  so  that 


(3.  32) 
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and  the  maximum  range  is 


R 

max 


The  corresponding  minimum  <j>  value  is  given  by  <J)^,  where. 


(3.  33) 


= tan 

T 


-1 


R, 


(3.34) 


Then  for  R^  » r^  we  have 


r 

t 


(3.  35) 


which  is  identical  to  This,  of  course,  follows  directly  from  the  geometry 

in  Fig.  5.  We  shall  use  this  model  in  the  next  subsection  to  include  the  effects 
of  shadowing. 


3.4  Shadowing  Effects 

The  evaluation  of  the  incident  power  density  upon  a given  aperture,  as 
developed  in  the  last  two  subsections,  did  not  take  into  account  the  possibility 
that  not  all  of  the  surface  is  irradiated.  When  the  surface  is  rough,  then  the 
very  nature  of  the  roughness  acts  as  a shadowing  mechanism  that  prevents 
the  incident  radiation  from  falling  on  a reflecting  surface.  Furthermore,  when 
the  radiation  is  scattered,  it  may  also  find  itself  shadowed  from  the  receiving 
point  as  a result  of  the  same  phenomenon. 
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Assume  that  we  have  some  arbitrary  surface  with  height,  ^(x,  y),  which 
is  a random  field.  Now  consider  the  geometry  shown  in  Fig,  6.  A wave  inci- 
dent at  angle  at  point  x^,  y^  is  scattered  at  angle  cj)^.  If  this  is  a specular 
type  of  scattering,  then  the  ray  in  the  cj)^  direction  is  the  specular  reflection. 
The  normal  n is  the  normal  to  the  surface  (x,  y).  The  height  of  the  surface  is 

given  by  ^(x  , y ) at  the  point  of  scattering.  Now  we  want  to  ask  the  question, 
s s 

what  is  the  probability  that  the  surface  y^)  is  irradiated?  What  chance 

is  there  that  the  ray  may  be  blocked  by  some  large  surface  perturbation? 

The  same  question  can  also  be  asked  about  the  scattered  ray.  To  answer 
this  question  we  shall  follow  the  work  of  Wagner. 

Wagner  has  derived  a function,  S(4).,cj)  ),  which  is  the  probability  that 
a ray  incident  from  angle  does  not  get  blocked  and  generate  a scattered 
ray  in  direction  cj)^,  which  is  also  not  blocked.  This  function,  called  the 
shadowing  function,  is  given  by 


= [1  - exp(-2(B.  + B^))] 

• [erf(vi>f  erf(v3)]/(4(B.  + B^)) 


where 


(3.  36) 


\ = [exp(-  erfc(vj^^/(4  ) 


(3.  37) 


V = I cot  ^ 


(3.  38) 
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Fig.  6.  General  geometric  scheme  for  shadowing. 
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where  k = i or  s and 


erf(x) 


exp(  -X 


2)  dx 


(3.  39) 


erfc(x)  = 1 - erf(x) 


(3.  40) 


Now  note  that  for  (J>.,  close  to  tt/ 2,  we  have  that  cot  (j)^  and  cot  c()g  are  close 
to  zero.  Thus,  and  Vg  are  close  to  zero.  Therefore, 


B 


k 


(3.41) 


so  that 


1 - exp 


Zy/n  V.  ZJff  V 
1 s 


Zjn  V.  V 
^ 1 s 

V.  + V 

1 s 


(3.  42) 


which  is  itself  quite  small.  Thus,  at  low  angles,  the  shadowing  is  quite  great. 

Furthermore,  even  if  just  one  of  the  cj)(s)  is  small,  the  corresponding  will 

be  quite  large  so  that  S(cj).,cj)  ) will  be  dominated  by  that  term  as  expected. 

X s 
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As  we  shall  note  in  the  next  subsection,  these  results  will  be  important 

in  modifying  the  channel  spread  function,  particularly  at  angles  away  from 

the  specular  points  where  <|).  or  <j)  approach  tt/2.  • Furthermore,  this  theory 

1 s 

is  also  useful  for  modifying  the  ideas  in  specular  multipath  theory. 

A more  recent  study  of  this  effect  has  been  made  by  Lynch  and  Wagner 
who  considered  both  multiple  scattering  and  energy  conservation  in  the  shadow- 
ing process.  Their  resulting  functions  show  a 3-dB  difference  at  lower  angles 
(1^)  and  would  probably  be  more  appropriate  at  grazing  angles.  Their  work 
is  based  upon  that  of  Wagner  and  also  Sancer. 

3.  5 Example  of  Channel  Spread  Functions 

In  the  previous  sections,  we  developed  the  theory  to  describe  the  sta- 
tistical characteristics  of  the  scattered  field  as  viewed  at  a receiving  aperture 
from  both  a flat  surface  and  a spherical  earth.  Extensions  to  other  arbitrary 
surfaces  are  possible  and  follow  as  a direct  application  of  the  theory.  In  this 
section,  we  evaluate  the  spread  function  for  several  different  geometries  and 
comment  on  the  limitations  of  the  theory  as  imposed  by  the  approximation 
used. 

The  spread  function  K(9,  <j>)  was  evaluated  for  a spherical  earth  as 
given  in  Equation  3.  25.  In  Figs.  7 through  11,  we  have  plotted  the  function 


lOlogio 


K(9,<j>) 

(PQ/4TTd^) 


where  d is  the  free  space  distance  between  transmitter  and  receiver.  All 
plots  are  for  sources  at  a height  of  1 km  above  the  earth’s  surface,  and 
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Fig.  8.  Channel  spread  (R  = 100  km,  s = 0.  5,  no  shadowing). 
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Fig.  9.  Channel  spread  (R  = 10  km,  s = 2.  0,  no  shadowing). 
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receivers  at  a height  of  10  m.  For  all  the  curves,  we  have  assumed  vertical 
polarization  for  the  transmitter  and  receiver. 

Figures  7 through  9 are  for  the  case  of  no  shadowing.  In  Fig.  7, 
the  spread  was  evaluated  for  s = 0.  5 and  a range  of  10  km.  It  can  be  noticed 
that  most  of  the  energy  is  concentrated  close  to  the  horizon  (low  elevation 
angles)  and  symmetrically  about  the  boresight  between  transmitter  and  re- 
ceiver. As  the  range  of  the  transmitter  (source)  increases  to  100  km,  the 
profiles  narrow.  There  is  also  some  redistribution  of  the  energy  density  not 
directly  noticeable  in  these  contours.  As  the  slope  increases  to  s = 2.  0 in 
Fig.  9,  a surprising  effect  occurs.  The  concentration  of  the  energy  in  the 
central  area  disappears,  and  sidelobe  behavior  arises.  This  effect  is  due 
to  the  increased  radius  of  curvature  term  in  the  scattering  cross  section. 

In  Appendix  I,  the  average  radius  of  curvature  is  given  by 


1 


1 


k^|q|^  D^(n^  • 


(3.43) 


However,  the  integral,  I^,  used  in  the  method  of  stationary  phase  was 


(3.44) 
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This  integral  was  obtained  by  assuming  that  the  distance  to  the  receiver 


from  the  scattering  surface  was  much  larger  than  the  dimensions  of  the  scat- 
tering surface.  Thus  using  Equation  3,  44  in  3.  43,  we  find 


and  the  receiver.  Therefore,  as  < R2  > ^ gets  large  and  becomes  com- 
parable to  this  distance,  the  basic  assumption  begins  to  break  down  and  the 
analysis  will  no  longer  hold.  Therefore,  for  very  large  s values,  the  value 
of  the  average  radius  of  curvature  can  become  excessive  for  certain 
values  of  7,  the  angle  between  the  z axis  and  the  normal  at  the  specular 
point.  Thus  it  may  be  appropriate  to  actually  limit  the  value  of  7. 


(3.45) 


Recall  that  is 


s 


4 


4 


(3.46) 


h 


Now  for  small  7,  those  close  to  the  line-of-sight  line  between  source  and 
receiver,  sec  7 is  almost  unity.  Then 


^4 

« 0.  1378  TT  ^ 


(3.  47) 


h 
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If  the  average  slope,  s,  is  unity,  then  for  large  correlation  lengths, 

cr  ^ can  be  comparable  to  the  distance  squared  between  source  and  receiver, 

clearly  a violation  of  the  assumptions. 

When  shadowing  is  introduced,  the  profiles  show  a drastic  change. 

Their  basic  shape  is  the  same,  but  there  is  a 10-dB  loss  in  amplitude  due  to 

shadowing.  This  effect  is  shown  in  Figs,  10  and  11. 

The  effects  of  shadowing  can  best  be  seen  by  evaluating  the  ratios 

of  power  densities  received  as  a function  of  range.  Specifically,  the  power 

2 

density  (watts/m  ) due  to  diffuse  multipath  is  given  by 


•TT 


K(0,  4))  dOdtJ) 


-TT 


2 

while  that  due  to  the  free  space  propagation  is  Pq/4tt  d . Thus,  in  Figs.  12 
and  13  we  have  plotted  the  ratio, 


as  a function  of  range.  Here  we  have  assumed  an  isotropic  receiving  antenna. 
This  is  in  effect  the  ratio  of  diffuse  to  direct  power.  In  Fig.  12,  we  have 
plotted  this  as  a f\inction  of  range  for  s = 2.  0 and  s = 0.  5.  Note  that  the 
rougher  surface,  s = 2,  0,  has  a higher  power  level  by  about  1 dB,  There 
is  also  some  increase  as  the  range  decreases.  The  difference  between  the 
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.DIFFUSE  TO  DIRECT  POWER  RATIO  (dB) 
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Fig.  13.  Ratio  of  diffuse  to  direct  power;  shadowing. 
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two  is  accounted  for  by  the  sidelobe  nature  of  the  s = 2*  0 case.  In  Fig.  13, 
we  have  plotted  the  ratio  for  the  case  where  shadowing  has  been  included. 

Now  note  the  drastic  difference.  First,  the  levels  are  10  to  20  dB  lower. 
Second,  there  is  a significant  change  with  range  as  expected,  i.  e.  , at  far 
ranges  the  incident  rays  are  close  to  the  surface  and  have  high  probabilities 
of  being  blocked.  Third,  the  curves  for  s = 0.  5 and  s = 2 have  been  reversed. 
Now  s = 0.  5 is  generally  6 dB  higher  than  that  for  s = 2.  0.  This  is  a res\ilt 
of  the  shadowing  occurring  in  the  sidelobe  region.  In  Fig.  14,  we  have 
separately  considered  vs  range  for  a slope  of  0,  3.  This  is  a rough  sea 
surface  value  as  discussed  by  Barton,  Again  note  the  significant  difference 
between  the  nonshadowed  and  shadowed  case.  Also  note  the  peak  in  the  non- 
shadowed  case  at  20-km  range. 

These  curves  clearly  indicate  that  for  these  surface  slopes  and  sensor 
geometries  the  diffuse  multipath  will  in  general  be  below  the  direct  path  level 
and,  as  will  be  shown  in  the  next  section,  have  little  effect  on  position  esti- 
mation accuracy.  Moreover,  when  other  surface  slopes  are  considered,  the 
levels  of  diffuse  to  direct  power  levels  are  about  the  same.  For  example,  in 
Fig.  15  we  have  plotted  vs  the  rms  surface  slope  s(s  = 0.  1 - 1.0)  for 

Zj  = 10  m,  = 1 km,  and  R = 30  km.  These  curves  are  for  both  shadowing 
and  nonshadowing  surfaces.  Clearly  for  low  values  of  s(s  = 0.  1),  the  diffuse 
power  is  greatest.  This  result  is  in  agreement  with  Barton* s recent  work. 

For  the  case  of  no  shadowing,  the  difference  between  s = 0.  1 and  s = 1 is 
5.  2 dB, 
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DIFFUSE  TO  DIRECT  POWER  RATIO  (dB) 
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Fig.  14.  Ratio  of  diffuse  to  direct  power 
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DIFFUSE  TO  DIRECT  POWER  RATIO  (dB) 


\ 


RMS  SLOPE  (S) 

Fig.  15.  vs  rms  surface  slope. 
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Section  4 


ANTENNA  PERFORMANCE  EVALUATION 

Diffuse  multipath  acts  as  a source  of  interference  to  angle  estimation 
processors  by  presenting  to  them  false  sources  of  energy  from  directions 
other  than  the  true  one.  Different  processors  have  different  performance 
characteristics.  In  this  section  we  consider  various  forms  of  processors 
that  result  from  different  models  for  the  received  signals.  Basically  there 
are  three  models  that  we  consider:  Gaussian  signal  in  Gaussian  noise,  an 
unknown  but  nonrandom  signal,  and  a known  signal.  Each  signal  model  results 
in  a different  processor.  We  also  consider  suboptimum  schemes  such  as 
monopulse  and  compare  it  to  the  optimum  schemes. 

The  purpose  of  this  section  is  to  consider  how  diffuse  multipath  affects 
the  ability  of  an  array  to  determine  the  angular  location  of  a source  of  radiation. 
By  considering  a wide  variety  of  processors  from  the  optimum  minimum 
mean  square  error  estimator  to  monopulse^it  is  possible  to  see  what  one 
gives  up  for  processor  simplicity.  We  evaluate  both  processor  form  and 
performance  for  a wide  variety  of  channel  models. 

In  this  section  we  will  consider  eight  different  processors  that  are 
used  to  estimate  the  angle  of  arrival.  Each  of  these  processors  results 
from  a different  model  for  the  signal  r,  that  results  from  the  time  processing 
of  the  r-f  outputs  of  the  N sensors  of  the  receiving  antenna.  Specifically  we 
assume  that  the  received  signal  at  r-f  has  the  form 
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r(t)  =yjz  Re[r(t)  exp{j  Zuf^t)] 


(4.  1) 


where  is  the  carrier  frequency  and_r(t)  is  the  N x 1 vector  of  complex 
envelopes  of  the  outputs  (in  phase  and  quadrature)  of  the  N sensors.  The 
N X 1 vector  r is  given  by 

(t)  dt  (4.2) 

which  represents  matched  filter  detection  of  r-f  burst  signals.  Other  more 
complex  time  variations  are  possible  and  have  been  discussed  elsewhere 
(see  Van  Trees  [2]).  The  N x 1 vector  r can  be  decomposed  into  three 
possible  components  as  follows 


f=#-b  m(|jL)  + + w 


(4.3) 


where 


b m(|jL)  = the  signal  component  due  to  the  direct  and 
specular  path.  |jl  is  the  angular  ^cation  of 
the  source  of  the  radiation,  and  b is  a scalar 
complex  envelope 

nj^  = the  diffuse  signal  component.  This  is  an 

N X 1 complex  Gaussian  zero  mean  random 
vector 

w = the  sensor  noise.  This  also  is  an  N x 1 

complex  zero  mean  Gaussian  random  vector 

and  E is  the  signal  energy, 
s 
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This  model  for  the  received  signal  is  a canonical  one.  The  variations 
in  processors  result  from  the  varying  assumptions  that  are  made  on  the  nature 
of  b.  Specifically  b may  be  a complex  Gaussian  random  variable,  an  unknown 
but  nonrandom  constant  or  a known  deterministic  constant.  Furthermore, 
processors  may  be  structured  with  and  without  n^^  being  present.  When  nj^ 
is  absent  the  result  is  a nonspread  processor. 

The  covariance  matrices  for  and  w are 


5„(^)  = E[n„nf] 


(4.4) 


K = E[w  w^] 


w 


(4.5) 


In  both  cases  they  are  Hermitian.  We  further  assume  that  (see  McCarty 

[4])  K ^ can  be  factored  as 
L J'  w 


= w^w 

W 


(4.6) 


where  W is  an  N x N complex  matrix  called  a whitening  matrix. 

To  determine  the  optimum  processor  for  each  possible  variation  we 
first  need  to  determine  the  likelihood  function.  This  is  determined  from  the' 
binary  hypothesis  problem  where  under  and  we  have 


H 


V 1 - ^ + w 


(4.7) 


^0=  £ = w 


(4.8) 
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Note  that  the  dummy  ox  hypothesis  contains  only  noise.  This  results  from 
the  assumption  that  n^^  arises  as  a direct  result  of  the  presence  of  the  spec- 
ular signal.  Thus,  when  that  portion  is  absent,  so  too  is  the  diffuse  component. 
Note  also  that  K^(|ji)  is  a function  of  p,  the  source  direction.  ^ 

The  likelihood  ratio  A(r,  ji)  is  given  by 

p(r|H  ) 

A(r,Hi)=  V...  . (4.9) 

p<£l  Ho' 

The  maximum  likelihood  processor  for  p then  is  the  one  that  uses  J)(r  ,[i) 
and  finds  that  value  of  \x  to  yield  the  maximum.  The  log  likelihood  function 
is  similarly  defined  as 

t(£,  Hi)  = in  A(£,hi)  . (4.  10) 

2 

Finally  for  unbiased  estimates  we  have  for  the  Cramer-Rao  bound  oner  , 


> 

-1 

E 

9M  J 

(4. 11) 


In  the  following  subsections  we  shall  use  this  canonical  model  and 
provide  statistics  on  b so  that  f(r,  and  the  C-R  bound  can  be  evaluated. 
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4.  1 The  Gauss-in-Gauss  Problem 


The  first  case  assumes  that  b is  a zero  mean  complex  Gaussian 
random  variable  with 

<r^=E[£rX*]  . (4,12) 


The  likelihood  function  for  this  processor  has  been  derived  by  McGarty 
[1  - 5],  and  is 


A(7  “ 

I^T 

exof-l  r^  K"^  R [R  + R rl 

1 1 
K 

1 V/  ' 

exp|^  r ^ W -J 

where 


(4. 13) 


“V’)  + R^(l*)  ,4.14) 

It  has  been  shown  elsewhere  that  + K^]  ^ r is  the  minimum  mean 

square  error  (MMSE)  estimate  of  the  N x 1 vector  bm(ii)  + nj^.  In  Fig.  16, 
we  have  used  this  fact  to  sketch  the  maximum  likelihood  estimator  in  the 
estimator-correlator  form.  The  processor  works  in  the  following  fashion. 

The  vector  r is  received  and  then  whitened  through  W.  Then  for  each  value 
of  |JL,  an  estimate  of  the  direct  signal  plus  diffuse  multipath  signal  is  made  and 
cross  correlated  with  the  whitened  signal  (denoted  as  As  |ji  is  scanned 
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Fig.  16.  Gauss-in-Gauss;  optimum  spread. 
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through  all  values,  f(£,  jj.)  is  recorded  and  that  value  of  |j.  maximizing  f(£,  IJ.) 
is  the  estimate.  In  general  for  practical  applications,  fixed  beams  are  formed 
and  suboptimum  processing  follows. 

The  performance  of  this  processor  could  in  theory  be  evaluated.  No 
known  exact  calculations  exist  due  to  the  dependence  of  on  p.  Parallels 

with  continuous  time  results  can  be  made  (see  Van  Trees  [Z],  Chpt.  7). 

In  a similar  fashion  we  can  consider  the  case  of  an  unspread  channel. 
For  this  case  we  set  £ identically  to  zero  (letting  = 0).  The  resulting 
likelihood  function  has  been  evaluated  by  McGarty  [1]  and  is 

2 

A(£,  p)  exp 


“Z"  £ 


K m(HL 
w 


- 1 

mV|i)  K 


r 

w — 


(4.  16) 


whe  re 


E 0-^ 
'T 


^ 1 A,, 

1 + E ov  m (|x)  K rn(|i) 
s b — w — 


(4. 17) 


This  form  can  be  put  into  the  estimator-correlator  form  by  noting  that  the 

A 

MMSE  estimate  of  the  scalar  b,  b,  is  given  by 


b=  cr^m^(p)  K r 


1 

w — 


(4.  18) 


62 


Thus  A(r,|jL)  becomes 

w’*’wm(M.)]  . (4.19) 

°b 

This  processor  is  diagrammed  in  Fig,  17,  Note  now  that  the  covariance  does 
not  depend  on  p.  Thus,  the  performance  bound  is  readily  determined.  Speci- 
fically 


= E 


2 2 2 
s “"T  “"b 


1 

s T 


9 

9p* 


B(p; 0)  + 


9p 


B(p,p')| 


p'  = o-i 


(4.20) 


where 


B(p,p')  = m (p)  K ^ m(p  + p') 

W 


(4.21) 
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18-6-16194 


Fig.  17.  Gauss-in- Gauss;  optimum  unspread. 
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algorithm  (see  McCarty  [4]).  This  model  is  useful  for  describing  inter- 


from  the  CR  inequality. 

In  the  preceding  processor,  it  was  seen  that  the  dependence  of 
on  p was  not  present  and  the  result  was  a simplified  processor.  A suboptimum 
processor  can  be  generated  by  assuming  that  the  specular  and  diffuse  com- 
ponents were  scalar  random  variables.  The  resulting  processor  would  be 
one  that  would  be  used  in  the  case  of  an  unspread  channel  where  merely  a 
power  or  energy  increase  was  observed.  Specifically  if  we  let 


where  I is  the  N x N identity  matrix  and  the  noise  spectral  height,  then 
the  suboptimum  processor  would  be 


ference  environments  in  the  ATC  context.  The  bound  oner  follows  directly 


(4.22) 


N 


0 


(4.  23) 


where 


r 


1 


(4.  24) 


I. 


1 


2 -,2  „ 2 ^ „ 
cr„  = N E cr,  + E 
B s b s 


K(p,  p')  B(p  - p' ) dp 


(4.25) 
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Careful  comparison  with  the  preceding  processor  indicates  that  this  sub- 
optimum processor  is  identical  to  that  of  the  nonspread  Gauss-in- Gauss 
problem  discussed  above.  Thus,  the  processor  is  diagrammed  exactly  as 
in  Fig.  17  with  the  appropriate  changes.  The  performance  of  this  processor 
is  given  by  the  bound 


■*; 

-1 

2 / 

1 

Al 

c 

B(m-) 

9m- 

K(m.')-^  B(hl')  d^.' 
1 9M' 

o 

(1 

- 

- 

(4.26) 


which  is  valid  only  for  azimuth  estimation.  Also 


1 


1 + 


K(m.')  B(m')  d^l' 


C = 


^l 


Nor 


1 + 


2 

O'. 


f-1 


1 2 

c 

n 
N 


K(^l')  dfi'  + 


(4.  27) 


and  K(|j.)  is  given  by  the  transformation 

' TT  / 2 


G(4,)  k(0,4,)  dc), 


K(n) 


\/l  - 


ir/i  /»2tt 


G(cl))  K(0,  (}))  d0  dcf) 


/-  TT  / 2 Jo 


(4.28) 


LS  a 


where  G((j>)  is  the  antenna  gain  in  elevation.  The  function  1 /yTT7i. 
Jacobian  resulting  from  the  transformation,  p = cos  0,  The  function  K(p)  is 
normalized  to  have  unit  area. 
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The  remaining  function  is  the  array  ambiguity  function  given  by  the 


expression 


B(tx)  = Im(p')  m1’(p'  + p)| 


(4.29) 


2 2 . . . 
wither^  the  interference-to- signal  ratio,  the  noise-to- signal  ratio  per  dipole, 

Kq  the  wave  number  (2tt/Xq),  and  d.  the  distance  from  the  array  phase  center 

to  the  ith  sensor  (dipole),  and  t denoting  the  complex  conjugate  transpose. 

2 . 2 

Furthermore,  we  can  relate  the  azimuth  variance,  to  a by 

y u. 


2 

^0 


2 

tan  0 


1 + 


2 cot  0 
sin  0 


2 

(T 

M. 


(4.  30) 


These  are  the  three  classes  of  processors  used  for  the  Gauss-in-Gauss 
problem.  They  all  entail  whitening  arguments  and,  except  for  the  first,  the 
performance  can  be  obtained  directly. 

A special  case  of  interest  in  the  Gauss-in-Gauss  problem  is  the  un- 
spread channel  with 


K = or  I 
w n 


(4. 31) 


and  a linear  array  with  spacing  d^  between  sensors.  For  this  case  it  is  easily 
shown  that 


(4.  32) 
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This  represents  a performance  benchmark  for  the  Gauss-in-Gauss  problem. 


4.  Z Unknown  but  Nonrandom 

A second  model  for  the  specular  source,  b,  is  as  an  unknown  but 
nonrandom  complex  variable,  i.  e.  , we  let 

b = (4.  33) 

where  A is  the  amplitude  and  (p  the  phase.  Since  both  are  unknown,  the 
likelihood  estimate  of  both  A and  <p  are  obtained  and  used  in  the  likelihood 
ratio.  Using  this,  the  likelihood  ratio  is  obtained  using  the  generalized 
likelihood  ratio.  Specifically, 


P(£l  H,,  A,  (p) 

A(£,  p)  = max  . (4.34) 

A,  <p  p(£|  Hq) 

This  technique  has  been  used  by  Hofstetter  and  DeLong  in  their  study  of 
amplitude  comparison  monopulse.  We  define 

K , = K + K (4.  35) 

1 n w ' 

K-  = K . (4. 36) 

0 w 
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Then  N beams  are  formed  by  taking  the  received  vector  r and  projecting  it 
onto  vectors  The  beam  gains  G.(p  - p')  for  the  ith  beam  are  defined  as 

G.(fJi  - p')  =£^(p')  m(p)  (4.37) 

where  p'  represents  the  known  pointing  directions  of  the  ith  beam.  We  further 
require  that  the  ) are  orthogonal,  i.  e.  , 


IN  i=j 

0 i/j 


(4.  38) 


We  shall  demonstrate  shortly  how  this  can  be  interpreted  in  the  sense  of 
sum-difference  monopulse.  Using  these  beams  we  obtain  M complex  measure- 
ments r.  where 
1 


r . = Ae^*^  G.(u)  + n.+w  . i=l,...,M 

1 gi  gi 


(4.39) 


where 


~t  - 


n . = g.  n 

gi  .^1  - 


(4.40) 


w 


(4.41) 
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Using  the  orthogonality  of^.,  it  can  be  shown  that  for  independent  w,  the  cor- 
responding 

= G w (4.  4Z) 

where 


w = 


g 


w 


gl 


w 


gM_ 


G = 


'V 

1 


gt 

--&M- 


Then  define 


K = G K,  G 
gl  1 


~ ~ ~t 

K ^ = G G' 
gO  0 


(4.43) 


(4. 44) 


(4.45) 


Using  this  we  can  obtain  the  log  likelihood  ratio  f(r,  jx)  as 


(4.46) 
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where  r is  the  M x 1 vector  of  beam  projections  and 


or 


S m(iJL) 


1 = G m(iJ.) 


Following  Hofstetter  and  Delong,  we  obtain  for  A 


and 


i 


- arg 


Using  these  in  the  likelihood  ratio  we  obtain 


(4.47) 


(4.  48) 


(4.  49) 


(4.50) 


(4.51) 
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Fig.  18  performs  as  other  processors  by  scanning  through  all  \x  directions. 

The  performance  of  processors  of  this  form  is  quite  involved  and 
have  been  discussed  at  length  by  Hofstetter  and  DeLong. 

4.  3 Known  Signal 

When  the  complex  number  b is  known,  a different  set  of  processors 
occurs.  Let  us  again  consider  the  case  where  r is  the  N x 1 sensor  outputs. 
Define 


lo  - 


(4.52) 


and  Kj  and  as  before.  Then 


?(r  , ji)  = In 


(4.53) 


Recalling  that  the  MMSE  estimate  of  n can  be  given  by 


(4.54) 
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18-6-16195 


Fig.  18.  Unknown  but  nonrandom. 
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Then?(r,  p.)  can  be  obtained 


(4.55) 


Using  the  whitening  decompositions  of  K and  defining 


±0-'^  = ■"lo 


(4.  56) 


the  resulting  processor  is  shown  in  Fig.  19.  The  first  portion  is  for  colored 
noise  estimations,  and  the  second  portion  is  for  known  signal  correlation. 

The  evaluation  of  the  performance  follows  directly  from  other  cases, 
however,  since  ^ is  a complex  function  of  p via  the  matrix,  an  exact 
formulation  becomes  quite  cumbersome. 

A special  case  worth  considering  for  the  known  signal  case  is  for 


K 

n 


= 0 


(4.57) 


K =0-^1 
w n 


(4.58) 
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Fig.  19.  Known  signal;  spread. 
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For  this  case,  ^(r,|j,)  becomes 


This  is  the  classic  beam  steered  case  found  in  radar  and  sonar.  To  evaluate 
the  C-R  bound,  we  first  note  that 


E 


Hi,  ji) 


(4.60) 


where  E = b b , and  d,  is  the  position  of  the  kth  sensor.  Let  us  further 

S tC 

assume  that  we  have  a linear  array  with  spacing  d^  between  elements.  Let 
arbitrary  starting  point  on  the  array  so  that 


<^k  = 


kd.^  + d 


M 


(4.61) 
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Substituting  into  the  expectation  sum  in  Equation  4.  60  we  find 


2 

k=l 


N/--\2  /oj\2  N 

k=l 


+(r)  + + V 


(4.62) 


Now  is  a free  parameter  that  we  can  use  to  tighten  the  C-R  bound.  We 
want  to  choose  dj^  so  that  the  sum  is  minimized.  This  yields 


cr^  > . 

12 

■-5^. 

O' 

n 

1 N(N^  - 1) 

(4.63) 


This  should  be  compared  to  the  bound  for  the  unspread  case  of  the  Gauss-in- 
Gauss  problem  that  had  the  additional  term  which,  for  high  vanished. 

4.4  Monopulse 

Monopulse  is  a special  case  of  known  signal  or  unknown  but  nonrandom 
signal  with  only  two  outputs  being  provided.  To  evaluate  the  performance  of 
monopulse  we  first  consider  the  unspread  channel.  The  two  outputs,  r^^  and 
r^,  represent  difference  and  sum  beams  respectively.  Let  and 

be  two  orthogonal  N x 1 vectors.  Then  define 
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s m((i)  +g^(li')  n 


(4.  64) 


and 


= ^ m(Hi)  + n 


(4.  65) 


Further 


l^(fi')  = m(^') 


(4.  66) 


^2^^  ) ^ (M-*  ) 


(4.  67) 


where  G^  and  G^  are  beam  steering  or  distribution  networks.  To  ensure 
even  and  odd  beams,  we  assume  N to  be  an  even  number.  For  a uniformly 
illuminated  array,  G^  is  ^ and 


^A  = 


■1 

0 


0 

■1 


0 ...  0 


1 . . . . 

+1. 


. •+! 


N/2 
N/2  + 1 


(4. 68) 


N 
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The  estimate  of  jlis 


M-  = 


1 

k 


(4.  69) 


whe  re 


k = 


9(M-  - M-')  ^(m-) 


(4.70) 


The  constant  k is  called  the  monopulse  slope.  For  a uniformly  illuminated 
linear  array  with  spacing  d^,  we  find  the  slope  to  be 


2tt  d 


k = 


0 N 
4 


(4.  71) 


For  values  of  |jl  - |jl'  <<  1 (small  off-boresight  positions,  |jl'  is  the  boresight 
angle)  we  have  following  Sharenson: 


2 

16 

^|ji,  mono 

/ 2ir  d ' 

2 

E 0-^ 

1 

•0 

1 

s b 

1 

1 X - 

1 N 

— Z~ 
% 

(4.72) 
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The  first  thing  to  note  is  that 


(4.  73) 

|i,  mono  |i,  optimum 

i,  e.  , the  monopulse  error  is  greater  than  that  for  the  optimum  processor. 
For  large  N,  the  difference  is 


2 

a 

|JL,  mono 

"2 

'^|i,  optimum 


1.  333.  . . 


(4.  74) 


In  general  this  is  insignificant,  but  it  does  show  that  monopulse  is  not 
asymptotically  optimum. 

The  final  case  of  interest  is  to  evaluate  monopulse  performance  for  a 
spread  channel.  If  we  define  the  difference  and  sum  gain  as 

- M-')  = (M-')  m(M-)  (4.75) 


G^|i  - |i')  = ^^p')  m(p) 


(4.76) 


and  the  noiseless  sum  and  difference  signals  as 


S = b G^P  - p') 

(4.77) 

A = b G^(p  - p') 

(4.78) 
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then  we  have  the  following 


= A + 


- |Jt")  dij."  + w 


(4.79) 


= S 


- |j.")  S'(|j.")  dfjL"  + i’2(M-')  w 


(4.  80) 


Defining  the  two  noise  terms  on  and  r^as  n^  and  n^, respectively,  we 
have 


= A + 


n 


(4.81) 


/» 

r 


S 


2 + 


(4.82) 


If  we  denote  and  n^^  as  the  in  phase  (with  A and  2)  components  of  these 

two  noises,  then  we  can  use  Sharenson!s  results  directly  to  prove  that 


(4.83) 


where 


] 


(4. 84) 
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(4.85) 


2 

'^Sc 


= E[nlj 


P = 


(4.86) 


Furthermore,  for  ||j.  - |jl' | <<  1 we  have 


Si  ^ 

s 


(4.87) 


Lsl  = 


k(|i.  - |j.*) 


(4. 88) 


By  straightforward  calculations  we  can  show  that 


..1 


p = Re 


K(|x  - pl")  G^(pl'  - |jl")  - |jl")  d|jL" 


(4. 89) 


°"ac  " / K(|j.  -|i.")  |g^(|j.'  - p")|  ^ d|j."  + Na 


'-1 


(4.  90) 


Sc 


K(pl  - pl")|g^pl'  - |jl")|^  dpi"  + Nol 


(4.91) 
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Using  these  in  the  value  for  cr  we  obtain 


Generally  for  |fi  - | <<  1,  we  can  neglect  the  last  term  and  also  the  next  to 

the  last  term.  Also 


(4.93) 


Using  these  approximations  we  have  for  the  case  of  small  channel  spread, 
with  respect  to  the  beam  width 


2 

0* 


cr  . + 

p, noise 


dp* 


(4.94) 
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Finally  by  normalizing  K(|j.')  by  the  diffuse-to-direct  signal  ratio 


(4.95) 


and  defining  the  rms  diffuse  multipath  spread  as 

1 


d^J.' 


da  1 


I 


K(p')  <V 


we  obtain 


(4.  96) 


2^2  2 
V ~ V»  noise  ^ diff 


(4.97) 


with 


2 _ 2 
^|j. , diff  ^D. 


(4.98) 


This  expression  corresponds  to  the  expression  given  by  Barton  and  Ward 

(P.  151),  where  we  have  incltided  the  Gg^.,  average  sidelobe  ratio,  already 
2 

in  K(|jl).  Thus  0^^  is  evaluated  directly  from  K(|i)  and  not  approximately  from 

2 

the  glistening  surface  approach.  The  reader  should  note  that  was  obtained 
in  the  last  chapter  and  was  shown  to  be  quite  small  (-10  dB  or  less)  for  most 
Casses  of  interest. 
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Also  one  should  note  that  0, 

da 

K(fi)  is  spread  over  the  beam,  i.  e.  , 


can  be  generalized  to  the  cases  where 


I. 


(4.99) 


K(p')  dfi 


2 


noise 


The  essential  factor  here  is  that  monopulse  performs  more  poorly  than  the 
optimum  processor;  yet  its  total  effect  is  small  compared  to  specular  errors 
noted  in  McGarty  [1], 

In  Table  4.  1,  we  summarize  all  the  processors  and  their  corresponding 
performance  values. 

4.  5 Performance  Comparisons 

The  analysis  of  the  previous  subsections  determined  the  structure  and 
performance  of  various  processors  that  may  be  used  as  angle  estimators.  In 
this  subsection  we  shall  concentrate  on  two  processors:  the  Gauss-in-Gauss 

suboptimum  processor,  and  the  monopulse  processor.  The  suboptimum 
processor  provides  a com^enient  lower  bound  to  array  performance.  It  is  a 
processor  that  scans  over  all  angles  and  utilizes  all  angular  information,  and 
in  so  doing  can  use  the  diffuse  field  and  extract  positions  information  from  it. 
The  monopulse  processor,  however,  is  a single-hit,  two-beam  estimator  for 
which  it  has  already  been  shown  that  diffuse  multipath  will  always  degrade 
performance. 
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TABLE  4.1.  PROCESSOR  SUMMARY. 
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TABLE  4.  1.  PROCESSOR  SUMMARY  (cont.) 
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The  suboptimum  processor  performance  is  given  by  Equation  4.  26, 
and  the  monopulse  performance  by  Equation  4.  92.  It  should  be  noted  that  both 
variances  depend  uponintegrals  ofthe  channel  spreadfunctionwith  the  array  ambi- 
guity pattern  on  the  difference  pattern.  These  integrations  can  be  quite  difficult, 
however,  by  introducing  Fourier  transforms  they  reduce  to  readily  evaluated 
finite  sums.  This  is  done  in  Appendix  IV,  where,  (f),  the  Fourier  transform 
of  K(|Jt)  is  introduced.  A specific  example  is  shown  in  Figs.  20  and  21  where 
K(^)  for  s = 0.  5,  R = 10  km,  = 10  m,  and  = 1km  are  shown.  In  Fig.  21, 
we  have  plotted  its  transform  Note  for  this  example,  the  transform  is 

effectively  zero  for  values  of  f greater  than  10.  In  Appendix  IV,  we  show  that 
the  necessary  performance  integrals  are  determined  by  evaluating  9^(f)  at  only 
those  points  determined  by  the  sensor  array  geometry.  This  means  that  in 
general,  (f)  need  only  be  evaluated  at  discrete  sets  of  points. 

The  behavior  of  the  suboptimums  and  monopulse  processors  vary 
greatly.  The  general  conclusion  regarding  the  monopulse  processor  is  that 
diffuse  multipath  always  degraded  performance.  This,  however,  is  not  the 
case  for  the  suboptimum  processor.  In  particular  by  looking  closely  at 
Equation  4.  26,  we  see  that  if  the  integral  of  the  CSF  is  positive,  then  it  is 
possible  that  the  performance  will  be  poorer  with  diffuse  multipath.  Yet  if 
the  integral  is  negative,  then  the  diffuse  multipath  has  information  about  the 
source  location  that  will  improve  location  accuracy.  To  see  this  in  some 
detail,  consider  the  following  example. 


H<A  (4.100) 

elsewhere 


Example:  Let  K(/Li)  be  given  by 


1 

Ta 


K(p)  = i 


I 0 
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DIRECTION  COSINE,  /i 


Fig.  20.  Channel  spread  function,  K(|jl). 
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t(f  )-TRANSFORM  OF  K(/i) 


FREQUENCY,  f 

Fig.  21.  Transform  of  K(|Ji);  X(f)* 
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Then  the  transform  is  given  by 


_ sin  2it  £A 

Iv  {A 


(4. 101) 


and  using  the  results  from  Appendix  IV  for  a linear  array  we  obtain 


(N-  |kl) 


(4.  102) 


where  we  have  used  the  symmetry  to  simplify  the  sum.  Canceling  and  letting 
X = kd/X  , we  define  a function  G(x)  as  follows 

G(x)  = (^(n  “ -^1 X sin(27T  X A)  . (4.  103) 

N 

Then  the  summation  can  be  written  as  -2  ^ G(x,  ) where  x,  equals  x at  the 

k=l  ^ ^ 

given  k values.  Now  the  sum  in  Equation  4.  102  is  evaluated  at  the  values  of 

X,  which  are  the  spacings  of  the  dipoles  normalized  by  the  wavelength.  In 

Fig.  22,  we  sketched  G(x)  for  two  different  values  of  A.  In  (a),  A is  very 

small  so  that  K(p)  is  quite  narrow  and  its  transform  is  very  broad.  Thus, 

G(x)  is  broad  and  all  the  samples  are  positive  so  that  Equation  4.  102  is  neg- 
2 

ative.  Thus,  cr  can  be  reduced  for  such  a channel  by  using  the  information  in 
the  diffuse  energy  because  it  is  coming  from  the  same  azimuth.  However,  as 
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G(X) 


-4-16413 


(a) 

FOR  A<<  l/N 


SAMPLE  POINTS 


FOR  A > 1/2  N 
Fig.  22.  Plot  of  G(x). 
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A increases  we  still  sample  at  the  same  x values,  but  G(x)  now  oscillates  so 

2 

that  Equation  4.  102  becomes  positive,  thus  increasing  This  can  also  be 

interpreted  as  the  diffuse  field  coherently  decoupling  the  array. 

This  example  can  also  be  applied  to  the  monopulse  processor  with  the 
result  that  performance  will  always  be  degraded.  The  basic  purpose  of  the 
example  was  to  show  that  the  suboptimum  processor  does  perform  significantly 
different  from  monopulse.  The  cost,  however,  is  that  the  suboptimum  processor 
must  scan  all  angles  looking  for  the  one  to  maximize  the  likelihood  ratio. 

In  Figs.  23  to  25,  (r^  has  been  plotted  for  boresight  pointing  vs  range 

for  s = 0.  3,  = 10  m,  = 1 km,  d/\  = 0.  5,  and  an  SNR  of  0 dB  per  element 

for  both  the  suboptimum  and  monopulse  processors.  An  isotropic  vertical 

aperture  gain  was  assumed  for  the  first  set  of  curves  (N  =1).  The  o*  ^ 

values  have  been  plotted  for  three  different  array  lengths  N = 20,  30,  and  40 

elements.  For  each  case  we  have  plotted  the  noise  limited  values  of  o*  . Then 

0*^,  with  diffuse  multipath,  was  plotted  for  a nonshadowed  surface  and  a shadowed 
0 

surface.  For  all  cases,  o*^,  for  the  optimum  processor  with  diffuse  multipath 
present,  is  less  than  the  noise  limited  value.  This  results  from  the  fact  that 
the  optimum  processor  has  used  the  information  in  the  diffuse  field.  Note  also 
the  dip  in  at  20  km  coincides  with  the  peaking  in  P in  Fig.  14.  The  results 
for  monopulse  are  also  plotted.  Note  that  they  are  consistently  higher  than  the 
noise  limited  value.  This  is  a result  of  the  fact  that  they  reflect  the  diffuse  field 
as  an  increase  in  the  background  noise  level. 

In  Figs.  23-25,  we  have  also  plotted  performance  vs  range  for  the  case 
where  vertical  aperture  is  being  used.  This  has  been  done  for  both  the  optimum 
and  monopulse  processor.  The  performance  curves  for  the  optimum  processor 
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AZIMUTH  STANDARD  DEVIATION  (deg) 
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Optimum,  Unshadowed 
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Fig.  23.  Array  performance  vs  range  (S  = 0.  3,  Z1  = 10  m, 

Z2  = 1 km,  N = 20). 
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'AZIMUTH  STANDARD  DEVIATION  (deg) 


Fig.  24.  Array  performance  vs  range  (S  = 0.  3,  Z1  = 
Z2  = 1 km,  N = 30). 


0 m, 
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-AZIMUTH  STANDARD  DEVIATION  (deg) 
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Fig.  25.  Array  performance  vs  range  (S  = 0.  3,  Z1  = 10  m, 
Z2  = 1 km,  N = 40). 


160 
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for  both  shadowed  and  unshadowed  surfaces  are  equal  to  the  noise  limited  value. 
The  monopulse  performance  with  vertical  aperture  and  a shadowed  surface  equals 
the  monopulse  noise  limited  value.  The  vertical  aperture  contains  twenty  ele- 
ments at  ha  If -wave  length  spacings  with  a sector  beam  formed  having  a cutoff 
beginning  at  5^  in  elevation.  The  vertical  pattern  was  formed  using  sin  x/x 
beams.  From  the  figures  we  can  see  how  greatly  the  addition  of  vertical  cut- 
off reduces  the  azimuth  error.  This  is  particularly  true  with  the  monopulse 
processor  with  unshadowed  surfaces  whose  performance  is  now  down  to  the 
noise  limited  value. 

A final  example  is  shown  in  Fig.  26  where  we  have  plotted  0*^  vs  s, 
the  rms  surface  slope.  As  s increases,  cr*  increases  for  the  suboptimum 
process  because  decreases.  The  opposite  would  be  true  for  the  monopulse 
processor.  It  should  also  be  noted  that  shadowing  does  play  a significant  role 
in  performance. 

To  conclude,  what  we  have  done  in  this  section  was  to  compare  the 
performance  of  the  suboptimum  processor  with  that  of  the  monopulse  processor. 
As  a result  of  this  comparison,  the  following  observations  can  be  made. 

1.  The  suboptimum  processor  can  have  an  order  of 
magnitude  performance  improvement  over  mono- 
pulse. However,  such  improvement  requires 
significant  hardware  and  signal  processing. 

2.  Performance  calculations  are  most  easily 
obtained  by  using  the  transform  of  the  CSF. 

The  transform  is  required  at  only  N -f  1 points. 
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AZIMUTH  STANDARD  DEVIATION  (deg) 


RMS  SLOPE  (S) 

Fig.  26.  Array  performance  vs  rms  slope. 
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3. 


The  performance  measures  provide  analytical 
means  for  comparing  processors.  Tradeoffs  in 
sensor  locations,  weighting,  and  phasing  can  be 
made  directly. 

4.  Surface  roughness  can  degrade  performance  with 
the  most  serious  degradation  for  monopulse  being 
at  s « 0.  1. 

5.  Vertical  aperture  with  vertical  cutoff  significantly 
reduces  the  effects  of  diffuse  multipath  and  should 
be  used  where  necessary.  Even  without  vertical 
cutoff,  a 20-foot  array  at  L-band,  with  s = 0.  3 
experiences  a maximum  diffuse  error  of  only  0.  24^, 

Experimental  work  has  been  performed  at  Lincoln  Laboratory  as  part 
of  the  DABS  program  to  evaluate  the  performance  of  monopulse  of  an  off  bore- 
sight  mode.  An  experimental  facility  has  been  constructed  (called  DABSEF) 
to  perform  this  measurement.  A planar  array  of  dipoles  is  used  to  form  sum 
and  difference  beams.  The  diploes  are  at  about  half -wavelength  spacings 
at  L-band  (1030  mHz  and  1090  mHz)  with  16  in  the  vertical  direction  and  32  in 
the  horizontal  direction.  The  horizontal  spacing  is  8.  36  inches,  and  the  vertical 
is  6.  24  inches.  The  vertical  aperture  elevation  pattern  is  shown  in  Fig.  27, 
and  the  sum  and  difference  patterns  in  Figs.  28  and  29,  respectively.  Using 
this  antenna,  radial  flights  from  DABSEF  out  to  ranges  of  50  nmi  were  flown, 
and  on  the  inbound  radial  the  azimuth  estimates  were  recorded.  Also  position 
locations  of  azimuth,  using  a boresighted  TV  optical  system,  were  recorded 
and  used  as  the  true  values  of  azimuth.  The  difference  between  the  monopulse 
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lNTENSITY{dB) 


Fig.  27.  Elevation  patterns  of  DABSEF  antenna  through  azimuth  boresight. 
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INTENSITY  (dB) 


Fig.  Z8.  Azimuth  cut  through  free  space  sum  pattern  at  elevation  of  0.  0 radian. 
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Fig.  29.  Azimuth  cut  through  free  space  difference  pattern  at  elevation  of 
0.  0 radian. 
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azimuth  estimates  and  the  optical  estimates  is  the  monopulse  errors.  These 
radials  were  flown  under  varying  ground  conditions  (foliage,  barren,  snow, 
etc,)  with  similar  results.  Also  the  antenna  could  be  tilted  down  5^,  shifting 
the  elevation  pattern  down  by  that  amount.  The  result  of  this  tilting  was  an 
increased  illumination  of  the  ground. 

In  Figs,  30  to  3Z,  we  present  the  results  of  several  of  these  flights. 

They  are  all  at  an  azimuth  of  308,7^  from  DABSEF  over  the  airfield  at 
L,  G,  Hanscom  AFB,  The  aircraft  heights  are  above  sea  level  (ASL)  values. 

The  phase  center  of  the  receiving  antenna  was  332*  ASL,  The  major  reflecting 
surface  in  the  first  Fresnel  zone  was  Hanscom  Field  at  130*  ASL,  The  ground 
was  clear  of  snow,  and  the  visibility  was  in  excess  of  30  nmi.  We  have  plotted 
both  the  power  measured  in  the  sum  beam  and  the  azimuth  error  as  noted  above. 
In  Figs,  30  and  31,  the  antenna  was  tilted  down  so  that  the  effects  of  the  diffuse 
multipath  were  accentuated.  In  these  two  cases,  the  instantaneous  azimuth 
errors  can  be  seen  to  be  large  at  closer  ranges,  (5  - 10  nmi)  as  predicted  by 
the  theory  in  Fig,  25,  As  the  range  increases,  the  errors  decrease  to  the  noise 
limited  values.  This  is  a result  of  the  shadowing  as  well  as  spreading  of  the 
diffuse  energy.  In  Fig,  32,  however,  the  azimuth  errors  are  of  uniform 
distribution  over  the  ranges  presented.  In  this  case  the  antenna  was  not  tilted 
and  the  vertical  cutoff  attenuated  the  effect  of  diffuse  multipath  as  predicted  by 
theory  (see  Fig,  25),  A careful  examination  indicates  that  diffuse  errors  are 
a factor  of  two  to  three  times  the  noise  limited  case  as  predicted.  The  results 
presented  here  appear  consistently  on  the  64  different  radials  flown  during  the 
experimental  program.  There  thus  seems  to  be  a clear  correlation  between 
theory  and  experiment.  Furthermore,  the  diffuse  errors  seem  to  be  independent 
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AZIMUTH  ERROR  (deg)  - SUM  SIGNAL  (dBm  ) 
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AZIMUTH  ERROR(deg)  SUM  SIGNAL  (dBm  ) 


Fig.  31.  DABSEF  data  (2K  , 308.7°,  antenna  down). 
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AZIMUTH  ERROR  (deg)  SUM  SIGN AL  (d Bm  ) 
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Fig.  32.  DABSEF  data  (4K  , 308.  7 , antenna  normal). 
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from  sample  to  sample  indicating  that  the  diffuse  field  acts  as  white  noise  and 
comes  from  different  locations  at  different  ranges.  This  also  is  consistent 
with  the  theory  of  Section  3. 
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SECTION  5 


CONCLUSIONS 

Diffuse  multipath  has  been  characterized  in  detail,  and  its  effect  on  the 
pointing  accuracy  of  a planar  array  has  been  evaluated.  The  model  used  in 
this  evaluation  is  one  of  many  that  are  found  useful  in  this  area.  It  provides 
the  antenna  designer  with  a function  that  is  readily  used  to  evaluate  array 
performance  when  such  multipath  is  present.  However,  as  discussed  in  the 
introduction,  this  is  but  one  of  several  such  models. 

As  a result  of  this  study,  the  following  conclusions  and  results  were 
obtained 

1.  The  Kodis -Barrick  model  is  useful  at  L-band  and  higher 
frequencies  because  all  the  constraints  that  are  required  for  the 
specular  point  approach  are  met. 

2.  The  channel  spread  function  was  evaluated  for  a spherical  earth. 

It  was  shown  that  this  function  was  obtained  directly  from  the 
scattering  cross  section.  This  function  for  a spherical  earth 
extends  the  form  developed  for  a flat  earth  as  discussed  in  McGarty 
[2].  The  spherical  earth  includes  the  effect  of  horizon  cutoff. 

3.  Shadowing  has  been  included  in  the  model,  and  the  resulting  channel 
spread  function  shows  a significant  change.  The  shadowing  effects 
developed  are  based  on  quasi-geometric  optics  that  coincide  with 
the  scattering  cross  section  analysis. 
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4.  The  evaluation  of  the  channel  spread  function  shows  that  diffuse 
energy  is  concentrated,  as  expected,  along  narrow  strips.  However, 
over  very  rough  surfaces,  sidelobes  seem  to  appear. 

5.  Without  including  shadowing,  the  total  power  density  arising  from 
diffuse  multipath  is  on  the  order  of  that  from  the  direct  path  alone. 
When  shadowing  is  included,  this  ratio  changes  significantly,  i,  e,  , 
diffuse  power  is  from  4 to  24  dB  below  the  direct  power.  This 
clearly  indicates  that  diffuse  power  should,  in  the  cases  analyzed, 
have  a minor  affect  on  azimuth  pointing  accuracy. 

6.  Antenna  pointing  errors  that  result  from  diffuse  multipath  are  small 
compared  to  noise  induced  errors  when  large  vertical  apertures 
are  used  to  attenuate  ground  reflections.  The  analysis  of  pointing 
errors  extended  the  present  theory  to  include  more  complex  surface 
models. 

The  numerical  results  used  in  this  report  are  for  aircraft  at  low 
elevation  angles.  At  these  angles,  diffuse  energy  is  small.  However,  at 
higher  angles  (<j>  > 10^),  as  may  be  found  in  satellite  systems,  the  diffuse 
energy  can  increase  significantly  and  the  curves  must  be  recalculated.  Yet, 
the  theory  presented  in  this  report  extends  readily  to  those  cases. 
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APPENDIX  I 


THE  SCATTERING  CROSS  SECTION 


The  evaluation  of  the  scattering  cross  section  depends  upon  an  analysis 
of  the  scattering  of  an  electromagnetic  wave  from  some  surface  S as  shown  in 
Fig.  I-l.  First,  center  a coordinate  system  on  the  surface  S.  If  we  consider 
the  point  P,  as  in  the  figure,  to  be  the  observation  point  with  coordinate  r^, 
then  Silver  (p.  132)  shows  that  the  field  at  that  point,  due  to  reflections  from 
the  surface  S,  is  given  by 


E(r  ) = , 

— — p 4ttU)C 


[(K  . V)V  + k K] 


exp(-jkr ) 


dA 


(1.1) 


where  S is  the  surface,  K the  current  density  on  the  surface,  k the  wave 
number.  The  differential  area  dA  is  the  quantity  dx*  dy*,  where  x*,  y*  represents 
the  coordinates  of  the  scattering  point  r^.  We  have  also  assumed  a mono- 
chromatic wave  with  frequency  CJ  and  that  the  surface  has  a dielectric  constant 
c.  The  V operator  operates  on  the  x*,  y*  coordinates.  The  range  r is  the 
distance  between  the  scattering  point  and  the  observation  point  P,  i.  e,  , 


and  r has  the  coordinates  x*,  y*.  Thus  V operates  on  the  r portion  of  r. 

^“S  ^“S 

To  obtain  the  total  field  at  the  observation  point,  one  must  sum  up  all  the  fields 
coming  from  surfaces  S that  comprise  the  scattering  plane. 


no 


V 


Fig.  I.  1.  Scattering  surface. 


Ill 


Since  the  coordinate  system  is  centered  about  the  scattering  surface 
as  we  are,  in  general,  interested  in  the  far  field  distribution,  we  can  make 
the  standard  approximations  for  r.  Namely, 


r=r-r+r'r  -2r*r, 
— p — p — s — s — p — s 


Dividing  through  hy  r^  • r^  and  assuming  that  l^gl  « 1 we  obtain 


(1.3) 


(I.  4) 


What  we  are  assuming  is  that  S is  some  small  bounded  surface  and  the  scatter- 
ing points  on  this  surface  are  at  a distance  small  compared  to  the  distance  to 
the  observation  point.  Taking  the  square  root,  we  further  approximate  r by 


r ^ 


(1.5) 


For  l^gl  small  compared  to  l£pl>  vector  r^  is  almost  collinear  with 
r - r . Thus,  the  direction  of  observation  is  the  direction  determined  by  r 
from  the  point  of  scattering.  We  define  this  by  the  vector 


— s 


(1.6) 


Thus,  the  field  due  to  surface  S is  coming  from  direction  k . The  total  field 

— s 

at  r is  the  sum  of  fields  coming  from  all  directions  or  equivalently  from  all 
P 

surfaces  S. 
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Using  this  simplification  we  can  make  the  following  approximation. 


exp(-jkr) 


exp  1 

f-j| 

r 

-F 

1 +j  k • r \ 
) ' — s — s / 

1 ^ 

Thus  the  term  in  Equation  I.  1 can  be  written  as 
(K- 

•(-j  '"Ispl) 


exp 


= -(K-  k^)k^ 


r 

-P 


exp(j  • r^) 


Using  this  in  Equation  !•  1,  we  obtain  for  the  observed  field 

exp  (-j  k|r  I) 


--P’  " . ^1  I 

^ 4iTue  r 
-P 


K exp(j  kj  r^)  dA 


(1.7) 


(1.8) 


(1.9) 


The  value  of  K,  the  surface  current  density,  must  be  known  to  obtain  a solution 
to  this  equation.  In  general,  such  a solution  is  quite  difficult,  but  if  the  Kirchoff 
conditions  are  used  (see  Jackson,  p.  283),  Equation  1.9  can  be  readily 
simplified.  These  conditions  use  the  field  incident  on  S as  the  value  for  K. 

Now  since  we  assumed  that  the  reflecting  surface  has  an  infinite  con- 
ductivity, we  then  have  for  a boundary  condition 


n X H.  = K 
— — ^inc  — 


(I.  10) 
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where  n is  the  normal  at  the  surface,  and  H.  the  incident  magnetic  field  vec- 
— ’ — ^inc  ® 

tor.  If  we  further  assume  that  the  incident  field  is  a plane  wave  with  electric 
field  vector 


— inc  " ^0  i -i  ' L'>  (I.  1 1) 

where  e is  the  polarization  vector,  and  the  incident  direction  vector,  then  we 

can  find  H.  (r)  directly  as 
—me  — ’ 

^ (i  iii  * ^0  (I.  1 2) 

Using  this  in  Equation  I.  9.  we  obtain 


) = 


4itW  €i 


exp 


r 

-P 


2 2 

This  can  be  simplified  if  we  note  that  k = cj  ep,. 

We  now  want  to  evaluate  this  integral  at  each  of  its  stationary  points. 
That  is,  if  we  let 


q=^-ki  (1.14) 

then  at  the  point  r,  which  has  a normal  such  that  for  the  given  values  of  k 
and  k^,  the  angles  of  incidence  and  reflection  are  equal,  we  have  a stationary- 
phase  point.  If  we  let  n be  the  normal  vector  at  that  point  and  let  ^(x,  y)  be 
the  z value,  then  the  scattering  surface  is  defined  by 
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f(x,  y,  Z)  = z - ^(x,  y)  = 0 


(I.  15) 


and  n is  given  by 


n=-2_ 

- Ivf| 


Using  Equation  I.  15  we  obtain 


n = 


9 i 


— e re 
9 y — y — z 


(I.  16) 


(I.  IV) 


and  e , e , e are  the  unit  vectors.  At  the  specular  point,  we  know  that  the 
— oc  ^y  ^z 

vector  q points  in  the  direction  of  n so  that 

i = ili  -ils  " 

We  now  seek  to  prove  this  fact  using  a stationary  phase  analysis  on 
the  integral  appearing  in  Equation  I.  13.  It  is  easily  seen  that  if  e^  is  the 
unit  vector  in  the  z -direction  and  n the  normal  of  the  surface  at  the  point 
(x,  y),  then 


dA  n • e = dx  dy  (I.  19) 

— “~Z 

This  is  nothing  more  than  the  Jacobian  of  the  transformation.  Define  the 
vector  h as  the  vector  cross  product 


h = k.  X e 
— — 1 — 


(I.  20) 
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Then  using  Equations  I.  19  and  I.  20  and  the  definition  of  q in  Equation  I.  13, 
we  obtain 


(1.21) 


The  vector  r is  given  by 


(I.  22) 


where  ^(x,  y)  represents  the  height  of  the  surface.  Thus  if  we  expand  ^(x,  y) 
about  some  particular  point  (x^,y^),  then  we  have 


Is  =^Ix+yiy  + 


i 8 6 , 1 4.  8 5 / I 


4.  1 8 , ,2  . 8 , , , , 

+ z — z <*  - *,  > + F3TT7  (*-*,)  (y  - y, ) 

<7  X 


4.  1 8 , ,2 
■^z--z'y-y,>  tz 

9 y ^ 


(I.  23) 


For  convenience,  define  the  following  constants 
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A = 7-  ^(x,  y) 


(I.  24) 


B.  = 


£ d xd  y 


e(x,y) 


X£, 


(I.  25) 


C = ^e(x,y) 

^ d y 


X , y 
£’  ^£ 


(I.  26) 


Then  the  product  q r will  equal 


Is " % 14)^  (s  It)  ^ 


‘Iz  14  ^£  % ^^^£’  '^£  ^ 


9 e 


T" 


[a-£  (x  - )^  + (x  - x^  ) (y  - y^  ) + (y  - y^  )^j 


+ . . . 


(I.  27) 


Now  the  method  of  stationary  phase  states  that  the  integral  in  Equation  I.  21 
vanishes  for  all  x,  y values  except  those  in  and  arovmd  the  points  where  the 
coefficients  of  the  x and  y terms  in  Equation  I.  27  vanish.  At  the  point  at 
which  they  vanish,  the  integral  is  nonzero  and  can  be  evaluated  (see  Born  and 
Wolf).  Thus,  the  stationary  phase  points  are  those  x^  and  y^  which  are 
solutions  to 
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(I.  28) 


+ q. 


a e - 

d X 


= 0 


q + q 


a e - 

ay 


= 0 


(I.  29) 


Since  the  region  where  Equations  I.  28  and  I.  29  hold  true  is  small  for  large 
values  of  |q|  , we  evaluate  the  functions  inside  the  integral  at  the  points 

' Thus,  they  are  constant  and  independent  and  can  be  removed  from 
the  integral.  The  value  of  the  integral  in  Equation  I.  21  evaluated  at  the  fth 
stationary  phase  point  then  is 


S£.p*  = - 

^ 4Trk 


exp 


Z ^0 


(-J  ’'lipi) 


r 

-P 


F 

—e 


n • e 
— f — z 


exp(-jqz(|4^j  +14^^ 


(a^  (x  - )^  + 2B^  (y  - ) (X  - x^  ) 


+ C (y  - y )^))dxdy 


(I.  30) 


The  term  F is 
—( 


F = ( 
—0  \ 


k - k^)  h X n 

- ~6  / - -I 


( 


Where  n is  the  normal  at  the  point  x , y , Using  the  fact  that 

—(  e ( 


(1.31) 


a e = _ X 
a X q 


(I.  32) 
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(I.  33) 


9 ^ = - S: 

q. 


we  first  note  that  using  Equation  I,  17 


qe  +qe  +qe  q 

X— X y— V z— «;  — 

n = — = 

/ 2 ^ 2 . 2 

-V  /q  + q + q 

\/^x  z 


(I.  34) 


at  the  stationary  point.  Furthermore, 


exp  (-q^  14  y?  + q^  ))) 

= exp(j  q^x^  + j q y^  + j q^z^  ) 


= exp(j  q • r^) 


(I.  35) 


Now  define  the  integral  in  Equation  I,  30  as  so  that  using  Equation 
I.  35  we  obtain 


- J 


F;  exp  k|  r^l  + j q • r 

^ ^0  2,  • ITI  “ 


) = 

P 47rk  ^ — z 


(1.36) 


r 

-P 


The  integral  can  be  readily  evaluated  using  standard  techniques  to  yield 


= 


2it  j 


q 

^z  I 


V7 


(I.  37) 


where 


D,  = A.C,  - B, 
I t I I 


(I.  38) 
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Now  recall  that  by  definition 


= i-  = lil£f  • iz 


(I.  39) 


since  n^  is  the  unit  normal  and  at  the  specular  point  is  in  the  same  direction 
as  q.  Using  this  in  Equation  I.  36  we  finally  obtain 


^ F,  exp  (-jklr^l  tjq-  r,) 

^ ^ (2;  • lal  lipl 


(I.  40) 


The  reflection  coefficient  for  the  surface  is  given  by 


(1.41) 


and  has  been  evaluated  by  Mitzner.  Define  the  vector  q as 

q = q/k  (1. 42) 

Using  this,  the  total  received  field,  as  a result  of  all  stationary  points,  is 

readily  obtained  by  summing  over  f,  i.  e.  , the  field  at  r^  coming  from  direction 

k due  to  a sources  incident  on  S from  direction  k.  is 
— s —1 


E(r^)  = 


V ^0 

exp 

Z "z" 

1 

(=1  *' 

|q| 

r D. 


~fT 


(2^ 


(I.  43) 


Now  the  scattering  cross  section  for  incident  direction  k.  and  received  direc- 


— 1 


tion  k is  defined  as  (see  Kerr,  p.  33) 


— s 
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(I.  44) 


(T  (k.,  k ) = 4tt 
—1  — ^ 


where  R = |r|  and  E[  ] represents  the  expectation  operator.  Using  Equations 
1.43  and  1.41,  we  readily  obtain 


O'  (k.,  k ) = E 
— i'  — s 


n n 


> y tt  R.(k.,  k ) R (k.,  k ) 
^ Z-!  i’— s m— i’ — s 


?=1  m=l 


exp(j  q • (£^  - £^)) 

|q|^  D (n^  • e (n  • e 

' -2. ' I m ^-4  —z'  m — z ' 


(I.  45) 


We  now  assume  that  the  specular  points  are  independent, and  distribution  such 
that  q • is  uniform  on  [O,  2ir].  Then 


0-  (k.,  k ) = E 
i’  — s 


y fei-hs’ 

I 


(I.  46) 


The  term  in  the  denominator,  |q|  ^ D.(n  • e )^,  can  be  shown  to  equal 

£ — z 


f ^ {^) 


I 9 - (9^1  ) 

Z"  1 Z"  \ 9 X 9 y/ 


.9  X 9 y 


n 4 


(«)  * (H)' 


/I  + 


(I.  47) 
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where  we  have  used  Equation  I.  34,  and  all  the  derivatives  are  evaluated  at 


x^,  y^.  Following  through  on  the  calculations,  we  observe  that 


1 

uMliML 

Iql'^  • 

'a  a '"e . (d  h Y 

A 

= < 


> 


f 


(I.  48) 


where  < R^  ^ ^ the  product  of  the  principal  radii  of  curvature  at  point  f. 

We  can  now  return  to  Equation  1.46  and,  assuming  that  the  specular 
points  have  equivalent  statistics,  we  can  write  Equation  1.46  using  Equation 
I.  48  as 


0-  (ki,  ^ ) = TT  E [n]  E[<  R^  R^  >]  R^(k.,  ^ ) (I.  49) 

where  we  have  noted  that  R(k.,  k ) does  not  depend  on  f.  Define  the  average 

—1  “~S 

number  of  specular  reflectors  as  n^,  and  the  average  of  the  product  of  the 
radii  of  curvature  as  Then 

^ iis  ^ iis  ^ 

which  is  the  fundamental  result  of  this  analysis.  We  can  also  define  the  scat- 
tering cross  section  per  unit  area  if,  instead  of  using  n^,  the  average  number 
of  specular  points,  we  use  n^  the  average  number  of  specular  points  per  unit 
area,  i.  e.  , 
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or  (k.,  k ) = Trn  .:^R^(k.,  k ) 
—1^  — s A ^ '--1^  — s 


(1.51) 


If  we  define  k.  and  k as  (see  Fig.  1-2) 
—1  — s 


k. 

“1 


rcos  0.  sin  cb: 
1 

sin  0.  sin  d>. 

1 

»-cos  <j)^ 


"cos  0 

s 

k = k^  sin  0 
— s 0 s 

.cos  d) 
^s 


sin  d> 
^s 

sin  6 
^s 


and 


60  = 0 - 0. 
s 1 


then  we  obtain  for  the  scattering  cross  section 


(1.5  2) 


(1.53) 


(I.  54) 


1 V ^ sec  7 


exp 


R (k.,k  ) 
—1  — s 


(1.55) 


where 


2 ^ .N  T 

R (k.,  k ) = (k.  X e X (k.  - k )) 
i’  — s —1  — —1  — s 


rk"kT-i)^(kkT.i) 

y — S — S — / V"™® — S — / 


(k.  X e X (k.  - k )) 
—1  - -1  — s 


(1.56) 
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with 


k 


— s 


(1.57) 


k. 
— 1 


(I.  58) 
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APPENDIX  II 


GEOMETRIC  CONSTRAINTS 

Consider  the  scattering  geometry  shown  in  Fig. 4,  Here  we  define  two 

angles  a.  and  p,  which  correspond  to  latitude  and  longitude.  Let  ^ be  the 

angle  from  the  horizontal  plane  of  the  sensor,  and  0 be  the  angle  from  the  normal 

n^  as  shown.  These  correspond  to  elevation  and  azimuth,  respectively.  The 

scattering  surface  is  at  angles  a , p.  The  radius  of  the  earth  is  given  by  R^^. 

E 

We  first  note  that 


e = 


(II.  1 ) 


Also  from  the  geometry,  if  we  project  a line  normal  to  the  plane  of  3 through 
the  scattering  point,  then  tan  :|>  can  be  shown  to  equal 

Zj^  + “ cos  a) 

sin-a 

E 

Similarly,  we  can  show  by  observation  that 
(■^  + = 


(n.  2) 


(n.  3) 


or 


c^s  = 5-4> 


+ a 


(II.  4) 


We  now  want  to  evaluate  the  vectors  and  s^. 
origin  O to  the  receiver.  Then 


rj  = (^1  + Rj.) 


Let  r be  a vector  from  the 


(n.  5) 
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where  n is  the  unit  normal  in  the  z directions.  Let  pbe  the  vector  from  the 
— z — 

origin  to  the  scattering  surface.  Then 


p = R^fcos  6 sin  an  + sin  6 sin  an  + cos  a n ) 

n E ^ —X  ^ — y — z 


(II.  6) 


Then  s is 


= £i  - £ 

~ ” P sin  a ri^  - sin  p sin  a ~ ^E 

In  a similar  fashion,  we  can  find  which  is 


(IL  7) 


£2  “ (^2  ^ ^E^  [sin  a'  + cos  a'  n^] 


(n.  8) 


where  a'  is  shown  in  Fig.  4.  Thus 

l2  = £2  - P 

= - cos  6 sin  a n + [(z„  + R„)  sin  a'  - sin  6 sin  a Rt^]  n 
^ E — X 2 E'  ^ E-*  — y 

m 9 

+ [(z  + R ) cos  a'  - cos  a R^]  n 

2 E E-"  — z 

We  now  want  to  find  the  angle  that  s makes  with  the  normal  n on 
the  surface  and  the  angle  0^,  which  is  the  angle  between  Sj^  and  s^  projected 
onto  the  plane  of  n^.  To  obtain  these  values  we  find  it  convenient  to  rotate 
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the  unit  vectors  and  n^  around  to  the  scattering  area  through  angles 

a,  p.  Using  standard  rotation  matrices,  we  find  the  resulting  vectors,  mj^, 
m ^1  and  m ^ to  be 


m.  = R n.  (i  = 1.  2,  3) 


(II.  10) 


whe  re 


*COS  P 

cos 

a 

sin  p 

-cos 

P sin 

a 

-sin  p 

cos 

a 

cos  p 

sin 

P sin 

a 

. sin  a 

0 

cos 

a 

(n.  11) 


This  yields 


"COS  p sin  a 
sin  p sin  a 
cos  a 


(IL  12) 


m. 


sin  p"* 
cos  p 
0 


(n.  13) 


cos  p cos  a 
-sin  p cos  a 
sin  a 


(IL  14) 
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In  this  case,  is  now  equivalent  to  n^  so  that 


m,  = n 
— 1 — s 


(II.15) 


Now  using  the  dot  products,  we  obtain 


4>i  = 


COS 


(11,16) 


and 


(II.  17) 


But 


12  „2  .2 
= R„  sin  a 
' E 


+ (R^(l  - cos  O' ) + z 


(II.  18) 


and 


js^l  = cos  p sin  a R^  + ((^2  ^ ^ ^E^^ 

2 

+ ((z  + R ) cos  a'  - cos  a R„) 


(II.  19) 
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Yet,  we  already  have  (j>  so  that  all  we  require  is  (}>.,  Thus 

s X 

2 2 

' s_2  “ ^ ^ 

+ sin  p sin  a[{^2  ~ ^ 

+ cos  «[(^2  ^E^  (11.20) 

We  also  have 


m 


ll 


TE]1 


= sin  (}>  cos  0, 


(11.21) 


m 


• g 

2 —2  = sin  d).  cos  0_ 

1 1 — 2 


il2 


—3  ■ -1 

111  I 


= sin  d)  sin  0, 
^s  1 


ai.  22) 


(II.  23) 


(II.  24) 
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Thus  combining  these  we  obtain 


(11.25) 


(II.  26) 


where 

2 2 

• Sj  = - cos  p sin  ci-  cos  a R^,  + sin  p sin  a cos  a R^, 

+ sin  a(rj  + “ cos  a))  (11.27) 

2 

• [(r^  + Rj-)  sin  a'  - sin  p sin  a Rj^]  + sin  a 
• [(r^  + Rj,)  cos  a'  - cos  a R^,]  (II.  28) 

• Sj  = - sin  p cos  p sin  a R^,  - sin  p cos  p sin  or  R^,  (11.29) 

sin  or  Rj.]  (II.  30) 
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The  angle  between  k.  and  k then  in  the  plane  of  n is  6 

— 1.  — ^ s 

63  = TT  + - 01  . (II.  31) 

Thus  (j>.,  (|)  and  0 are  the  values  that  correspond  to  (j>.,  cj>  and  60  in  the  last 

X S S X s 

section. 

Having  c|)g  and  0^,  we  have  enough  information  to  describe 
which  we  know  depends  on  only  these  three  parameters.  Thus  we  make  the 
equality 


<!>,:  <l>i) 


(XI.  32) 


Also,  both  and  have  been  evaluated  in  Equations  II.  18  and  II.  19, 

respectively.  The  only  remaining  term  is  dA.  From  the  geometry  we  know  that 


dA  = R sin  a 


da  dp 


(II.  33) 


From  Section  3,  we  know  that 


dp  = d0  (U.  34) 

We  now  want  to  find  da  in  terms  of  d<f>.  Recall  that  Equation  II-2  gives  <J>  in 
terms  of  a so  that 
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d tan  (()  = sec  (}>  d(() 


= (esc  a - esc  a cot  a 


R, 


+ 1 da 


(II.  35) 


But; 


4.  2^  ^ , 
sec  <j>  = tan  <J>  + 1 


zi  + RjjCl  - cos  a ) 


2 


Rj,  sin  a 


+ 1 


z j + R^(l  - cos  a)^  + R^(l  - cos  a ) + R^  sin^a 

;2  ."2 

sin  O' 


(II.  36) 


Yet  from  Equation  II.  18,  we  see  that  the  denominator  is  nothing  more  than 
I s j I Therefore,  using  Equation  II,  35  in  Equation  IL  34,  we  obtain 


do'  = 


. 2" 

R^  sin  a 
E 


d({) 


2 / zj 

CSC  a “ CSC  a cot  — + 1] 


(II.  i7) 


Combining  with  Equation  II.  3 3,  we  obtain  for  dA 


dA 


£ii 


2 

CSC  a 


CSC  Q cot 


■ft  • ). 


d9  d<}) 
sin  a 


(II.  38) 


which  is  the  desired  result. 


133 


APPENDIX  III 


SHADOWING  EFFECTS 

Consider  the  geometry  shown  in  Fig.  III.  1,  Here  we  have  depicted  a 
ray  coming  from  infinity,  incident  at  x,  y equal  to  zero.  The  height  of  the 
surface  is  given  by  ?(x),  and  it  is  sufficient  to  consider  only  one  dimension 
variation.  Let  the  curve 

z = + (tan  il>)K  (III.  1) 

represent  the  incident  ray  where  is  the  angle  between  the  x,  y plane  and  the 
normal  (x  axis).  Note  that  (//  = tt/Z  - <J>^.  Furthermore, 

= ?(0)  . (III.  2) 

Now  define  ) as  the  conditional  probability  that  the  point  is 

illuminated  given  and  (the  slope  at  , y^)  with  the  incident  ray  defined 
as  in  Equation  III.  1.  Then  by  definition 

S(0Uo>  illuminated  I (HI.  3) 

= p[e  will  not  cross  + tan  ip  x 

for  any  X >0 

lim  s(<// Up,  T) 

T -*  00 
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Fig.  III.  1.  Geometry  of  shadowing  effects. 


135 


where 


S(^Uq,  t)  = P[^  will  not  cross  + tan  ip  x 


for  X e(0,  t)Uq  ] 


Now  define 


q(r)  dr  = P[^  crosses  the  ray  inx€(r,  T+dr) 


ii 


does  not  cross  in  x e (0,  t); 


^0^ 


Then 


^0’ 

= ^q5  [1  - q(T)  dr] 


or 


(III.  4) 


(III.  5) 


(III.  6) 


(III.  7) 
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This  then  implies  that 


s(i//||^.  t)  = C exp 


q(u)  du 


(III.  8) 


Now  as  T -►  0 we  have 


lim  S(i,||  |';t)=0 

T-n) 


if  >tan  i// 


(III.  9) 


or 


lim  t'I=  I if  ejj  <tan  <//  . (IK.  10) 


That  is,  if  the  slope  at  the  origin  is  greater  than  that  of  the  incoming  ray,  then 
the  point  will  be  shadowed.  Therefore, 


C = u_,  (tan  JJI 

where  u ) is  the  unit  step  function.  We  really  want  S(i/y),  the  probability 
that  any  point  will  be  irradiated  by  a ray  with  angle  (//.  This  is  the  uncondi- 
tional probability  of  Equation  IIL  8.  But  we  assume  that  we  know  the  joint 

probability  densities  of  ^(x)  and  ?'(x).  Let  this  be  p*  -,(u,  v).  Then 

s>  s 


S(0)  = 


(u,  v)  du  dv 


(lU.  12) 
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Now  Equation  III.  12  can  be  solved  if  we  had  q(r  ).  This  is  difficult  to 
in  general,  so  we  approximate  it  by 

q(r  ) dr  ~ P[^  crosses  + tan  ip  x in  xe(r,  t + dr  ) 

with  slope  > tan  Uq,  ^q] 

This  probability  has  been  evaluated  by  Wagner  and  is 


q(r  ) = 


; — 2—7: 

2^  o-Q  0-^ 


tan  Ip  = V. 


0 


exp 


(^0  + Po 


2.1  / 2.^ 


where 


O-Q  = 


obtain 
(ni.  13) 


(IIL  14) 

(ni.  1 5 ) 

(IIL  16) 
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From  our  discussion  in  Section  3,  we  know  that 


E[^(x)  ^(x')]  = (T^  exp^-rVo’^j  (III.  17) 

2 2 

where  here  r = (x  - x')  . Thus 

(m.  18) 


and 


2 


cr 


1 


(III.  19) 


Furthermore,  we  know  that 

E[e(x)^|(x)]  = 0 


(III.  20) 


for  this  given  correlation  function  which  implies  that  ^ and  are  independent. 
Now  Equation  in.  18  can  be  used  in  Equation  III.  8 and  then  in  Equation  III.  12 
to  obtain  S(ip).  However,  this  is  for  the  monostatic  problem.  For  the  bistatic 
problem  as  shown  in  Fig.  6,  we  have  a (//.  and  or  equivalently  cf).  and  cf)^. 
Then  the  probability  that  we  see  a wave  scattered  to  ^ ^ from  cj).  given  and 
becomes 


4^' -O’^o) 


(ni.  21) 
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However,  for  the  cases  of  interest,  (j).  and  represent  scattering  into  two 
different  and  independent  directions.  Thus 


(III.  22) 


Therefore,  following  the  procedure  outlined  before,  we  have 


<t>g)  - 


Pt  Pt  t(v)  du  dv 

^0  ^0 


(III.  23) 


This  integral  has  been  evaluated  by  Wagner  to  yield  the  result  of  Section  3. 
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APPENDIX  IV 


PERFORMANCE  BOUND  CALCULATIONS 

In  general,  the  integrals  appearing  in  Section  4 are  difficult  to  evaluate. 
However,  if  we  assume  that  K(fj,)  is  such  that 

K(p)-0,  p = + 1 (IV.  1) 

then  considerable  simplification  arises.  This  is,  in  general,  the  case  as 
discussed  in  the  last  section.  We  now  evaluate  the  integrals 


K(hl') 


B(p'  ) d(JL' 


(IV.2) 


K(hl'  ) B(p')  dp' 


(IV.  3) 


and  the  integrals  appearing  in  the  evaluation  of  the  performance  of  the  mono- 
pulse processor. 
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These  integrals  can  be  related  to  the  integrals  of  their  Fourier  transforms 
via  Parseval's  theorem.  Namely,  let 


-6  (f) 


exp(-jZTrfn)  B(n)  d^JL 


(IV.  4) 


7t(f) 


exp(- j Zirfp,)  K()jl)  dp. 


(IV.  5) 


Then  we  can  write  for  Equations  IV.  2 and  IV.  3 


K(p) 


B(n)  dp.  = 


(ZTTf)^;e(f)^(f)  df 


(IV. 6) 


K(p)  B(p)  dp  = 


;<</(f);^(f)  df 


(IV.7) 


Now  note  that 


m (p«)  m'  (p*  + p) 


N 

^ exp(-jZirpdj^A) 

i=  1 


(IV.8) 
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so  that 


N N 

B(jjl)  = EE  exp(j  2tt  |i(d.  - dj^)A) 

i=l  k=l 


The  Fourier  transform  of  this  is  obtained  quite  readily  as 


N N 

^ (f)  = ^ ^ 6(f  - (d.  - dj^)A) 

i=l  k=l 


Using  this  in  Equations  IV.  6 and  IV.  7 we  obtain 


B'(M-)  dfx 


and 


N N 


K(fx)  B(hi) 


i=  1 k=  1 


(IV.  9) 


(IV.  10) 


(IV.  11) 


(IV.  12) 


143 


Thus,  knowing  the  Fourier  transform  of  the  channel  spread  function 


at  the  specified  number  of  points,  the  integrals  in  Equations  IV.  6 and  IV,  7 are 
directly  evaluated. 

A special  case  of  interest  is  that  of  a uniformly  spaced  array.  Namely, 


k + 1 


k = 1,  2, 


N 


(IV.  13) 


For  this  case,  the  double  sum  in  Equation  IV.  10  can  be  shown  to  be  equivalent  to 


N 

(f)  = E (N  - |k|)  6(f  + k d/X)  . (IV.  14) 

k=-N 


The  corresponding  integrals  for  a uniformly  spaced  array  are 


(IV.  15) 


N 


K{^,)  B(hl)  dp  =V  (N  - |k|)^fe'j 


(IV.  16) 
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Thus  for  symmetric  we  must  evaluate  the  transform  at  only  N + 1 loca- 

tions. This  is  a considerable  savings  over  direct  evaluation  of  the  integral. 
One  can  also  note  that  the  evaluation  of  the  second  partial  derivative  of  B(M-) 
at  fJ.  = 0 is  equavalent  to  using  a K(p)  equal  to  6(|Jl).  Thus 


6(p') 


9p-' 


B(p')  dp' 


= - (2 


N N , ,2 


(IV.  17) 


i=  1 k=  1 


For  the  case  of  a uniformly  spaced  array  we  have 


9p 


B(p) 


N 


p=  0 


= - (2tt)^  ^ (N  - |k|) 

k=-N'  ' 


(IV.  18) 


Now  Equations  IV.  15,  IV.  16,  and  IV.  17  can  be  used  to  evaluate  tr 

M 
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For  the  case  of  monopulse  performance  with  a linear  uniformly  spaced 


array  we  have 


) = - 


k=l 


exp|j  k 2tt  (fj.  - fj.'  )j 


(IV.  19) 


Using  standard  substitutions,  we  can  show  that 


4/N 


G^(hl  - ) = 4 


sin  l-a ^ 2tt((JL  - (JL*  )J 


"X 

sin  1*2^  2tt(^l  “ |x*  )j 


(IV,  20) 


Following  the  use  of  transform  techniques  as  discussed  previously,  we  can 
now  directly  evaluate  the  integral  appearing  in  the  monopulse  bound.  Namely 


K(hl)  \ 


2 


dfJL 


(IV. 21) 
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where  is  the  transform  of  G^(fJi).  Using  the  previous  techniques 

readily  shown  that 


K(hi)  |g^(hi)|^  dHi 


N/2 


= N 'H{0)  + 2^  (N  - 3k)9i(^j 
k=l 


These  expressions  can  now  be  used  directly  to  evaluate  performance, 
that  in  all  cases  the  transform  need  only  be  evaluated  at  N + 1 points. 


> 


, it  is 


(IV.  22) 


Note 
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GLOSSARY  OF  TERMS 


A 

antenna  area 

b(0.  4>) 

complex  envelope  of  time  invariant  portion  of  diffuse  field 

b(0.  4>,t) 

complex  envelope  of  diffuse  field 

c 

velocity  of  propagation 

d. 

1 

distance  to  ith  dipole 

d 

o 

E[  ] 

interdipole  spacing 
expectation 

erf  ( ) 

error  function 

E 

energy  of  transmitted  signal 

E 

— sc 

scattered  field 

f(t) 

complex  envelope  of  transmitted  signal 

G(0.  <//) 

antenna  gain 

g^(0.4)) 

difference  antenna  gain 

G^0.  4>) 

sum  antenna  gain 

Bessel  function 

i(_) 

current  density 

k 

monopulse  slope 

k 

plane  wave  vector 

Z ' 

components  of  propagation  vector 

K(0) 

channel  spread  function;  azimuth  only 

K(0.  4>) 

- channel  spread  function 
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K^(r,  r') 

covariance  of  diffuse  signal 

m(|jL) 

- array  delay  vector 

N 

number  of  array  elements 

n 

normal  vector 

n(t) 

complex  envelope  of  receiver  noise 

"a 

average  number  of  scatterers 

N 

o 

~T 

noise  spectral  density 

P 

o 

power  of  transmitter 

R 

ground  range 

R{?) 

reflection  coefficient 

radius  of  earth 

range  to  horizon 

Rrp 

range  from  aircraft  to  antenna 

r 

range 

r 

point  in  space 

r(t) 

complex  envelope  of  received  signal 

^0 

scattering  plane  displacement 

^1 

location  of  receiver 

Lz 

location  of  transmitter 

s 

scattering  surface 

^(4*  j ' 4*2^ 

shadowing  functions 

s(t) 

transmitted  signal 

s(t) 

complex  envelope  of  transmitted  signal 

S(£) 

complex  envelope  of  time  invariant  diffuse  signal 
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s(r,  t) 

- 

complex  envelope  of  time  variance  diffuse  signal 

SNR 

- 

signal-to-noise  ratio 

t 

- 

time 

^0 

- 

signal  emission  time  of  source 

- 

transmission  time  from  transmitter  to  receiver 

t 

s 

- 

transit  time  of  scattered  path 

var  ( ) 

variance 

Xi) 

- 

antenna  phase  center  location 

^1 

- 

antenna  height 

^2  1 

^2/ 

- 

location  of  aircraft 

^z) 

z(X,  Y) 

- 

random  surface 

r(r,  r') 

- 

Greens  function 

7 

- 

noise-to-signal  ratio 

A 

- 

difference  signal 

V 

- 

del  operator  (gradient  or  divergence) 

6(k  - k') 

- 

delta  function 

e. 

1 

- 

dielectric  constant  of  ith  scattering  plane 

V 

- 

diffuse  multipath  to  free  space  signal  ratio 

/\ 

9 

- 

azimuth  estimate 

- 

beamwidth 

X 

- 

free  space  wavelength 

M 

- 

direction  cosine 

Mb 

- 

beamwidth 

151 


?(x.  y) 
? 

Pd 

2 

% 

""h 

2 

‘"I 

2 
cr . 

1 


n 


cj)^;  4>i) 

2 

2 

cr 

|J- 

4>i 

‘t’s 

5^,, 


0 


Ri 


CO. 


random  surface  height 
interference-to-signal  ratio 
diffuse  reflection  coefficient 
sum  signal 

channel  attenuation  coefficient 
standard  deviation  of  surface  height 

diffuse  signal  energy 

power  in  diffuse  field 

scattering  cross  section 
conductivity  of  ith  scattering  plane 
correlation  length  of  surface 

noise  covariance 

signal  energy 

scattering  cross  sections 

variance  of  azimuth  estimate 

variance  of  p 

elevation  (incident  or  ith  signal) 
elevation  of  scattered  signal 
phase  difference 
multipath  phase 
phase  due  to  range 

carrier  frequency 
array  ambiguity  function 
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